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hermal History of the Earth with consideration of the 
Variable Thermal Conductivity of its Mantle 


H. A. Lubimova 


(Received 1957 September 30) 


Summary 


Data about the abundance of the radioactive elements and the 
thermal conductivity of the Earth are in favour of the view that its 
initial temperature T(r) must have been below the melting tempera- 
ture. The results of the calculations of To(r) are presented for the 
Earth considered as being formed through accumulation of particles 
of a gas-dust protoplanetary cloud. 

The redistribution of radioactive elements is considered from 
uniform distribution to a layered one. 

The matter of the mantle is considered as a dielectric possessing 
not only molecular conductivity but a radiative transfer of heat in 
the high temperature range. 

The Green’s function method and hydraulic analogy method are 
used for the calculations of the solution of the thermal conductivity 
equation. 

It is shown that the increase of effective thermal conductivity due 
to a radiative transfer of energy at a great depth does not produce any 
essential losses of inner heat. The thermal conductivity of the Earth’s 
upper layers is largely determined by its molecular component decreas- 
ing as the temperature rises and thereby forbidding considerable losses 
of inner heat. That leads to the conclusion that the Earth’s thermal 
history may be regarded as secular heating of its interior with simul- 
taneous slight cooling of its outer layers. 


1. Introduction 
Several methods for the estimation of the temperature inside the Earth exist 
‘at present. The possibility of determining the temperature from the degree of the 
increase of electrical conductivity with depth and from the velocities of seismic 
"waves has been shown. However, only an investigation of the thermal history of 
H 115 











116 H. A, Lubimova 


the Earth can give an answer to the question whether the Earth is undergoing the 
process of cooling or heating. A number of problems of the Earth’s evolution are 
connected with this question. 

Experimental data about the heat flow lost to space, the Earth’s age, density 
distribution in it, content of radioactive elements in rocks, thermal conductivity 
and thermal capacity of rocks are used for the study of the thermal history of the 
Earth on the basis of the heat conduction equation. Moreover, definite suggestions 
regarding the initial distribution of temperature and the function of the distribu- 
tion of the heat sources in the Earth are also introduced. 


2. Initial conditions 


The temperature of fusion, or of the solidification of the Earth’s matter, is most 
frequently accepted as the initial temperature To(r). On this assumption however 
any heat generation in the deep interior of the Earth, whence no heat is escaping, 
should be excluded. Otherwise the Earth mantle’s temperature would now exceed 
the melting point and the mantle could never become solid independently. The 
heat transfer in the interior is assumed to take place mainly by conduction as it is 
in the paper of Jacobs & Allan (1954) or by some more extensive mechanism. 
The small thermal conductivity of the solid outer layer would forbid considerable 
losses of inner heat and heat would accumulate inside the Earth, raising the tem- 
perature higher than the melting point temperature. 

There is no reason to believe that the mantle is quite free from radioactivity 
at present. ‘The data about the heat flow on the oceanic and continental areas 
having granite basalt layers of different thickness (Revelle & Maxwell 1952) give a 
direct indication that such sources are present deep inside the Earth’s mantle. 
The increased amount of data concerning the content of U, Th and K in various 
rocks and meteorites leads to the conclusion that radioactive elements are spread 
all over the Earth (Birch 1956, Verhoogen 1956). Even in the iron meteorites, 
with the composition of which the core of the Earth is usually identified, a quite 
definite concentration (Davis 1950, Paneth, Arroll & Jacob 1942, Starik & 
Shatz 1956) of them was discovered. In addition to that, the amount of iron in 
the Earth’s core is often supposed now to be much less than was believed earlier 
(Bullen 1957). According to the hypothesis of Lodotchnikov (1939) and Ramsey 
(1949) the composition of the Earth’s core differs very little from that of silicates. 
Hence we cannot completely exclude radioactive elements even from the Earth’s 
core and cannot assume 7 o(r) = 7T(melting). 

Cosmogonical ideas represent one of the foundations for conclusions about the 
initial temperature and the distribution of radioactive elements. According to 
modern ideas the Earth and the planets were formed as a result of an accumula- 
tion of particles of a gas-dust protoplanetary cloud. The most complete explana- 
tion of the main features of the structure of the solar system from the point of view 
of this hypothesis is given in the papers by O. Y. Schmidt and his collaborators 
(Schmidt 1954, 1957, Levin 1956). These papers also lead to the conclusion that 
To(r) was less than 7(melting). The initial temperature of the Earth was caused 
by the following three processes: (a) impacts of particles against the surface of the 
growing Earth; (b) compression of the matter inside the Earth under the pressure 
of the growing outer layers ; (c) radioactive heat, which could accumulate during the 
process of the Earth’s formation if this process was not very short. 

Effect (a) was estimated by Safronov (1954, 1958) for two limiting schemes of 
the Earth’s growth: (1) at the expense of large particles, which are not mutually 
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colliding; and (2) at the expense of small inelastically colliding particles. If the 
variations in the particle velocity are not taken into aecount the maximum increase 
of temperature is about 200° K in the case of the first scheme and about 1200° K 
in the second scheme at a distance of about o-8R from the Earth’s centre. 

Taking into account the variation of the velocity of particles in the first case, 
which is most probable, the maximum temperature increase will only be 50—-100° K. 
The period of the formation of the greatest part of the Earth (about 97 per cent of 
its mass) is about 0-23.109 years (Safronov 1958). A correction for radioactive 
heating for this period is about 700° K (if the Earth’s age is about 5.10% years). 
The temperature rise due to the increased pressure in the process of the Earth’s 
growth can be estimated from the formula (Lubimova 1955a) 


— exp| —(— . ~)| (1) 
2a\ por _—p? 

where 79 is the temperature of the matter before its compression, a is the co- 

efficient characterizing the dependence of Gruneisen’s constant on the volume, 

and po is the density of the matter before its compression. If all the above effects 

are taken into account the maximum temperature rise in the process of the Earth’s 

formation is about 1200° K in its centre. 


T.(r)°K 


T melting 





Fic. 1.—Distribution of the initial temperature 7,(7) in the Earth formed as a 
result of accumulation of gas-dust particles of a protoplanetary cloud. 


Curves of the distribution of the initial temperature To(r) are given in Figure 1. 
The hatched and the dotted curves were calculated without taking the variations 
of the velocity of particles into account and neglecting the radiogenic heat in the 
period of the Earth’s growth (the hatched curve corresponds to scheme (1) of the 
Earth’s growth; the dotted—the scheme (2) ). These corrections were taken into 
account in the preparation of the continuous curve. As it is seen from Figure 1, 
To(r) is much less than the temperature of fusion (thick line), as given by Uffen 
(1952). Thus the Earth formed as a result of an accumulation of particles of the 
gas-dust protoplanetary cloud would not become fused. 

The distribution of the radioactive elements in the Earth was usually assumed, 
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when the thermal regime of the Earth was calculated, not to vary in the process of 
the Earth’s evolution (Slichter 1941, Urry 1949, Jacobs & Allan 1954). 

However, modern geological ideas lead to the conclusion about the permanent 
formation of the Earth crust, the gradual growth of the continents by means of a 
partial remelting of an initially uniform meteoritic mantle and carrying of the easily 
fusible sialic constituents from the inside to the surface of the Earth, as a result of 
which the distribution of the radioactive elements would have been changed 
(Rubey 1951, Jeffreys 1956, Wilson 1954, Magnitsky 1955). 

From this point of view the equality of the heat flow on the oceans and the 
continents may be understood. The total amount of sources must be the same in 
these parts of the Earth, but underneath the oceans the matter is less differentiated, 
the heat sources being distributed more uniformly than in the continental regions. 

The suggestion about the uniform distribution of the radioactive sources of heat 
during the initial stage of the Earth’s evolution follows naturally from the hypothe- 
sis of its origin as a result of an accumulation of particles from the protoplanetary 
cloud. In the process of the Earth’s evolution there had been a gradual passage 
from a uniform distribution of sources to the modern layered distribution. 

The period of the existence of the uniform distribution of sources can be 
assumed approximately equal to the difference between the values of the age of the 
Earth as a planet ((4-5).109 years) and of the age of the Earth’s crust (3.109 years). 


3. Heat conduction with constant conductivity 


3.1 Uniform distribution of sources 

A solution of the equation of heat conduction containing the function of sources 
H(r, t) at constant x, c, p was obtained by means of Green’s function using the 
method of reflections (Lubimova 1952, 1956a, 1956b) for its construction. This 
method, developed in the mathematical papers of Tychonoff (1937) and Tychonoff 
& Sumarsky (1951), leads to a rapid convergence of series. After some trans- 
formations an expression for the temperature inside a spherical body is obtained in 
the form of a difference of two items having a simple meaning. The first depends 
upon the time and represents the temperature resulting from the heat generated 
in a given volume. The second, in the form of a rapidly converging series, charac- 
terizes the influence on the temperature of heat escape towards the surface*. 
For the thermal diffusivity of rocks k = x/cp = 0-01 cm?/s («—thermal con- 
ductivity), and for time intervals less than 6.10° years it is sufficient to take the first 
term of the series, if the body radius is greater than 1600 km (for the Earth and the 
Moon in particular). From 4 to 16 terms are taken into account (Lubimova & Starkova 
1954) for the radii of 200, 100, 50 km (the dimensions of the largest asteroids). 

The temperature inside the Earth for a uniform distribution of the sources 
and of the zero initialt and the surface temperature is given by the formula 

* The obviously applied Laplace transformation method (Lowan 1933) or of the division of 
variables (Allan 1955) led to extremely slowly converging series. Their summing up has become 
possible only since the electronic computers were introduced (Jacobs & Allan 1954). It may be 
shown that the method of automodel solutions leads to the same slowly converging series. 

+ The influence of the initial temperature T(r) can be taken into account separately by means of 


an additional term representing the solution of the homogeneous equation of thermal conductivity 
for the cooling globe, deprived of internal sources of heat: 


R 
- | T(r’).G(r,t,r’,0)dr 
0 


where G(r, t, 7’, r) is the Green’s function (Lubimova 1952, 1956). 
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t t 
I R (R- 1) 
T(r, t) = — | H(r)dr — —— | H(r). O| —————|dr 
(1) af (7) = ©) oot (2) 
H(t) = >, Halts)exp(—Aat) @(x) = 1—erfx (3) 


where A, is the decay constant of the element a, and H,(to) is its rate of heat genera- 
tion fp years ago. From formula (2) it follows that the region of heat escape towards 
the surface, and the inner region where all the emitted heat is spent on heating may 
be distinguished in the temperature field of the Earth. In the internal region 
r < R—C(kt)** the second term is negligible, the temperature being deter- 
mined by the first term, increasing with the rise of t. 

In the surface layer (r < R—C(kt)*) the second term becomes comparable 
with the first one. Some portion of heat flows from this region towards the surface. 
The depth of the region of escape is proportional to (k#)*. For a usually 
accepted value of k = 0-01 cm2/s it equals about 1000 km and may reach the 
centre after 101° years only. It follows that the cooling of the Earth as a whole 
cannot begin earlier than after 1019 years. If k = 0-1 the region of heat escape 
will exceed 2g00km already after 2.109 years of the Earth’s formation 
(Lubimova 1953). Heat generation below the region of heat escape, as well as the 
changes of the physical properties of matter which are going on there, do not in- 
fluence the heat flow observed at the Earth’s surface. 

The obtained solution permits us to calculate the thermal flow 


oT 


= K*— 


or \r=R 


lost by the Earth. Results of a calculation of g for various values of the mean 
concentration of the radioactive sources and of the different values of age are given 
in Figure 3. The contents of U, Th and K in the Earth were taken equal to their 
mean abundances in meteorites according to the two sets of values A and B 
obtained earlier (Lubimova 1952). According to the value A one gramme of the 
Earth contains 0-93.10-8§ g of U, 0-67.10-1%g of AcU, .2-8.10-8§g Th and 
1-95.10~? g of K4°, According to the value B: 1.10~? g of U, 0-72.10~9 g of AcU, 
3.10~? g of Th, 1-95.10~? g of K49, A new value, C, is obtained from recent data 
(Starkova 1955): 5:2.10-8g U, 3-7.1071%g AcU, 21.10-8g Th, 0-9.10-’g of 
K40, This value C gives the curve of heat generation, which differs little from 
the curve of variant A for the last 2.109 years and practically coincides with it in 
the more remote past (Figure 2). All these values are within the intervals limiting 
on the one hand the minimum value introduced recently by Urey (1956) and on the 
other hand a maximum value, which may be deduced from the new data by Starik 
& Shatz (1956) (Figure 2). 

We see from Figure 3 that the surface heat flow has decreased during the last 
2000-3000 million years. However, conclusions about the cooling of the Earth as a 
whole, such as were supposed by Urry (1949), are not correct. The Earth’s interior 


* C is a coefficient depending upon = H,(tg). 
a 
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which is below the region of the heat escape is getting heated. Only the surface 
layers are cooling. 

The concentration of radioactive elements, the age of the Earth and the lost 
heat flow are mutually connected. The calculated values of the flow g are within 
the range of the measured values (1-2 + 0-6).10~® cal/cm?s, if concentration A and 
the Earth age of (4-5).109 years are accepted (see Figure 3). 


* 
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H(t): 107? cal/cm*® sec. 
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Fic. 2.—Heat generation, according to different values. 


Formulae obtained for the temperature assuming a uniform distribution of heat 
sources are applicable to the description of the thermal regime of bodies having 
much smaller masses than the Earth, as for example the Moon and the asteroids. 
Differentiation of the matter of these bodies may be neglected. It may be considered 
in some approximation that these formulae describe also the thermal regime of 
the Earth below the oceans, where there is a very low differentiation of matter. 

The features of the thermal field and of the heat flow of the Earth pointed out 
above are maintained, as it will be shown below, both for a non-uniform distribu- 
tion of heat sources, as well as for a variable coefficient of thermal conductivity. 

We suggested that the change from a uniform distribution of heat sources to 
their layer distribution took place 3.109 years ago. The temperature distribution 
established by the moment of the differentiation was accepted as the initial one. 
The total content of U, AcU, Th, K in the Earth corresponded to the variant A. 
The process of redistribution of sources is supposed to have taken place only in the 
mantle and the crust of the Earth; the region of the core was not involved. As a 
result of this the heat generation in the 20 km crust increased ten times (up to 
1-3.10713 cal/cm?s), while in the mantle up to a depth of 2gookm it has cor- 
respondingly diminished 8 per cent. 


3.2 Layered distribution of sources 


The solution of the equation of heat conductivity for a layered distribution of 
heat sources was also obtained by means of the Green function with the aid of the 
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method of reflections (Lubimova 1956). Curves calculated from these formulae 
are given in Figures 4 and 5. The changes that took place before the differentia- 
tion were not taken into account in the curves of Figure 4. Curves of Figure 5 
give the resulting temperature distribution*, in which the influence of the heat 
discharged in the stage of the Earth formation and the stage of uniform distribu- 
tion of radioactive sources were taken into account. From these curves it is seen 
that the high concentration of the sources in the Earth’s crust cause but an insigni- 
ficant increase of its temperature, because the largest part of the heat generated in 
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Fic. 3.—Surface flow of heat for the different Earth’s ages and contents 

of U, Th, K. Continuous line is the variant A, hatched line is the 

variant B. Straight lines a and b are the limits of the interval of observed 
values q. 


the crust was rapidly lost in space. The increase of the heat flow g was 20-30 per 
cent as compared with the flow at a uniform distribution of sources. At present 
the heat flow equals 1-26.10~® cal. per cm?/s which is close to the observed one. 
A comparatively small increase of q is explained by the fact that the dimension of 
the region of heat escape exceeds the thickness of the Earth’s crust many times. 
We see that about 2.109 years ago the heat escape to the surface began to exceed 
the total amount of heat generated in the crust and carried out from the deeper 
layers. A lowering of the temperature of the upper layers has begun as a result of this. 

The curves of Figure 5 confirm the conclusion that the thermal histories of the 
upper and internal layers of the Earth are different, which was arrived at from the 
simple instance of a uniform distribution of the heat sources. An accumulation of 
heat inside the Earth continues even at present. The same conclusion seems to be 
made recently by Jacobs & Allan (1954). ‘The temperature of the base of the mantle 
has increased about 200° in the last billion years. The temperature of the mantle 
in this place is about 4000°C at present. Thus the thermal history of the Earth 
may be represented as a secular heating of its interior, rather than its secular cooling. 

* For the calculation of the curves of Figure 5 the content of the radioactive elements was taken 


according to the variant A, the Earth’s age equal to 4.10° years and the initial temperature 7)(r) 
in accordance with the dotted line of Figure 1. 
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Fic. 4.—Temperature distribution near the surface of the layered Earth. 
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Fic. 5.—Temperature distribution in the Earth calculated at the assump- 

tion of constant coefficient of thermal diffusivity k m o-o1 cm?/s. Thin 

continuous line is the temperature at present, hatched line is the tempera- 
ture 2-10; years ago, thick line is the melting temperature. 
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Another property of the temperature field of the Earth is the cooling of the 
upper layers, which are inside the zone of the heat escape. 


4. Thermal conductivity of the Earth 


The above conclusions depend essentially upon the thermal conductivity of 
the Earth, which determines the rate of the losses of heat by the Earth. Let us 
discuss this question. 

According to its conducting properties the mantle of the Earth resembles 
dielectrics most closely. The thermal conductivity of dielectrics is determined by 
the vibrations of their crystal lattice. Therefore the estimates of the changes of 
thermal conductivity x, obtained according to the data about electric conductivity 
(McNish 1937) proceeding from the Wiedemann-Franz law, are applicable only 
in the region of the metallic core of the Earth. 

In order to obtain a better estimate of the thermal conductivity of the mantle 
one should proceed from theory of the lattice thermal conductivity of solids. 
Debye (1914) has shown that the main property of a crystal lattice leading to its 
finite thermal conductivity is the anharmonicity of vibrations of its atoms. The 
thermal motion of atoms must be represented as a totality of mutually interacting 
plane waves (phonons). 

Peierls (1929, 1955) investigated this question from the point of view of 
quantum mechanics, taking the discrete character of the crystal lattice into account. 
Only the collisions of the phonons of a discrete lattice lead to a finite thermal 
conductivity, as a result of which the summary wave vector is changing; this 
change is equal to a period of the inverse lattice. These are the so-called ‘‘Umk- 
lapp” processes. At high temperatures, T > 6 (@—the Debye temperature) the 
“Umklapp”’ processes are as frequent as the ordinary collisions. At low tempera- 
tures T' < @ the probability of such processes is small, the free path of the phonons 
being large; therefore the scattering of the phonons on the lattice defects, impurities, 
or mosaic blocks has a considerable influence on the thermal conductivity of 
dielectrics. In the present paper one high temperature case is of interest, when 
such effects are not appreciable. Debye and Peierls took the third-order terms into 
account in the expansion of the potential energy in a power series with respect to 
the displacement of atoms from their equilibrium positions. A similar dependence 
« « T-1 was obtained by them in the high temperature region T > 6. The 
same dependence was obtained according to semi-empiric estimates by Dugdall & 
MacDonald (1955) and Kontorova (1956), who also restricted themselves to the 
approximation of the third order. 

Pomeranchuk (1941, 1943) has shown that the temperature dependence of 
thermal conductivity obtained by Debye and Peierls takes place only for bodies 
having the special form of dispersion of velocity of thermal vibrations and special 
dependence of the velocity upon the direction. In the general case of an isotropic 
body anharmonicity alone can give no thermal resistance and the integral deter- 
mining the thermal conductivity* 

ws (Cyigp 
m1 © { puplp Ey (4) 
will be infinite unless fourth-order terms are included in the expansion of the 
potential energy, according to the powers of the displacement of the atoms from 


* The symbols in formulae (4) are: /,—the free path of a phonon, C,—the specific heat of a 
phonon, u,—its group velocity, dp = dp,dp,dp,; p/h—is its wave vector. 
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equilibrium positions, that is, collisions in which four phonons are participating 
are included. 

The peculiarities of such collisions had been examined by Pomeranchuk (1941, 
1943), who obtained a solution of the kinetic equation for this case. 

The following expression was obtained for xj: 





/ Mu2 \5/4 
~— kNua{ 


(5) 


where N is the number of elementary cells in 1 cm, k is the Boltzmann constant, 
M is the mass of a cell, a is the lattice constant, u is some mean sound velocity, 
which is determined from u, the velocity of longitudinal elastic waves and u, 
the velocity of transverse elastic waves by the formula 


I 2 ae 
ue 3 ( ru ug ) 

The dependence of the thermal conductivity on pressure is determined 
according to formula (5) by means of values of u and a. 

The matter of the Earth possesses isotropic velocities of sound; its temperature 
is higher than the Debye temperature. This gives ground for supposing that for- 
mula (5) may be applied to a description of the phonon component of thermal 
conductivity of the mantle. 

Apart from ordinary conduction by means of lattice vibrations, the radiative 
transfer of heat, due to the scattering of electromagnetic waves, is possible in dielec- 
trics at a very high temperature. 

It has been established experimentally (McQuarrie 1954) that the thermal con- 
ductivity of ceramic materials of the MgO, AlgO3, BeO types decreases as 1/7” 
where n = 1, for temperatures T < 1500-1800°K. However at T > 1500-—1800° 
it begins to rise rapidly. The same was established for glasses (Ginsburg 1930). 
At these temperatures along with the usual molecular thermal conductivity the 
radiative transfer is introduced. 

Since the mantle of the Earth consists of substances similar to ceramics and 
glasses, the role of the radiative transfer of heat inside the Earth must be taken into 
account (Clark 1956, Preston 1956). Indeed the investigations of the olivine 
spectra (Keller, Spots & Briggs 1952, Runcorn 1956) show that the maximum 
radiation at temperatures prevailing in the Earth will take place in the interval of 
respective transparency. 

As was shown by Van der Held (1952) the effect of radiative transfer in solids 
may be taken into account by introducing into the equation of thermal conduc- 
tivity, instead of the usual coefficient of thermal conductivity, the coefficient 


consisting of the sum x = «1+ «2 where, taking the radiative transfer of heat into 
account, «2 is 


16 n2oT3 
K=—° 


3 € 
« is the coefficient of the extinction of the material, m is the refractive index, a is 
the Stefan-Boltzmann constant. Such a method of account follows from the trans- 
formation of the usual equations of radiative transfer in a grey body, provided that 
the solids possess a considerable absorption coefficient. There are grounds for 





(6) 
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considering the matter of the Earth as a grey body, because of a high spreading of 
the olivine spectra. 

As a whole the thermal conductivity of the Earth’s mantle must be determined 
by the sum of phonon and radiative thermal conductivity: 


ueut = 16 n2a T3 
=. .° (7) 








The coefficient of proportionality B in (7) was determined by us from the con- 
dition that the molecular thermal conductivity x; must be equal to the normal ther- 
mal conductivity of rocks x; = 0-006 cal/g.s.deg C on the Earth’s surface at u 
= 36 km/s, p = 2-7 g/cm’; T = 277°K. As it is seen from formula (7) the first 
term is decreasing with the rise of temperature, while the second term is increasing. 
This circumstance is very essential for the study of the thermal history of the 
Earth. 

Clark found that with decreasing coefficient of extinction from ¢« = © to 
«= 100—107! cm the temperature of the Earth’s interior must fall 3 times. 
However the decrease of the molecular thermal conductivity «x; with temperature, 
which takes place in the upper layers and prevents considerable amounts of heat 
being carried out from the Earth’s interior (owing to which no rapid cooling of the 
Earth can take place (Lubimova 1957)), was neglected by him. 

The calculation by Clark was carried out on the assumption that temperature 
distribution is stationary in the Earth. In consequence of this he was obliged to take 
a 100 times less radioactive heat generation in the Earth, in order to satisfy the heat 
flow observed at the surface. 


5- Heat conduction with variable conductivity 


Let us consider the non-stationary equation of thermal conductivity 


PP: Ms [a+ maya ~|+ 2 H, exp (—Agf). (8) 


ot r2 Or 


It becomes nonlinear when the changes of the coefficient of thermal conductivity 
k = k1+«2 and of density p(r) are taken into account. The method of hydraulic 
analogies of Lukyanov (1939) was used for the solution of this equation with 
the aid of the hydro-integrator. This apparatus permits the simulation of the pro- 
cess of thermal conductivity for variables x, c, p in one, two and three-dimen- 
sional space. The change of the amount of heat in an elementary volume V;, cm, 
per time unit equal to dO/dt = cpV;,dT/dt, is simulated by changing the quantity 
of water dO/dr = w,dA,/dr where A), is the level of water in a volume V;,, a, 
the area of its cross-section in cm?. A similar analogy is established for the heat 
flow q = K(T;,—T;,_,) from volume V;, to V,_,; and the water flow g = r,4_; 
(A;,— A,_1) between V;, and V,_,. K and 7,,_, are the coefficients of exchange, 
which are the values inverse to the resistances : the thermal value 
I dea 

Revs K k.F 
(x«—thermal conductivity, F—section of the area crossed by the heat flow gq, 
d;. ;-y—the distance between the volume centres) and the hydraulic value: 





1/pk,k-1 = Tk, k-1- 
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If the method of dimensions is used, the criterion of similitude is obtained, from 


which the scale of time m, and the law of simulation of heat generation follows from 
the formulae 


ceRk, x-1 : cpmagmy 


A 
k,k-1 = 


T ’ 
WkPk, k-1 OK .m, 


where /;,,-, is the hydraulic resistance determining the water flow, which simu- 
lates the heat generation; m, and my are the scale of temperature and the scale of 
intensities of the internal sources respectively. 

The precision of calculations by means of the hydro-integrator depends upon 
the degree of the division of the model into blocks. In our calculations the error did 
not exceed 7 per cent. It may be reduced to 1 per cent. 

The calculation of the thermal history of the Earth, taking the variation of the 
thermal conductivity according to the depth into account, was carried out according 
to the scheme of its evolution given in the first part of the present paper. But the 
thickness of the Earth’s crust was accepted as 40 km. The mean heat generation 
in the Earth’s crust was calculated according to the data on the content of U, Th, 
K in the granites and basalts given by Birch (1954), which coincide approximately 
with those given by Bullard (1954). The ratio of granites and basalts in the crust 
was taken as 1:2. This leads to heat generation 3.1071 cal/cm’s at the present time. 
The initial temperature To(r) was taken to correspond with the continuous curve 
of Figure 1 deduced on the ground of recent and more precise estimates. The 
mean abundance of radioactive elements was also taken from more precise data 
according to the estimate C (see above). The age of the Earth was taken equal to 
4°5.10° years (Patterson, Brown, Tilton & Inghram 1953)*. The redistribution of 
radioactive sources of heat was supposed to take place about 3-0.10® years ago. 

According to the data of a comparative analysis of the density distribution in the 
Earth and other planets of the same chemical composition (Koslovskaya 1953, 
Levin 1956) the differentiation of the matter in the Earth affects only the upper 
part of the mantle. We supposed in this connection for the last model that the 
carrying out of the sources into the Earth’s crust took place in the upper part of the 
mantle to a depth of about 1400 km. The concentration of radioactive elements in 
the impoverished part of the mantle will correspond to a dunite one. Below 1400km, 
including the region of the core, the source content was supposed similar to what 
it had been before the differentiation, corresponding namely to the mean of their 
content, C, in the iron and stone meteorites. 

The phase transition on the boundary of the core must have been going on 
when the mass of the Earth reached 0-8 of its present mass (Ramsey 1949). We 
believe therefore that the core had metal conductivity already in the stage of uni- 
form distribution of sources. The thermal conductivity of the Earth’s core was 
calculated according to the law of Wiedemann-Franz. 


k\2 
k= 3(=) cf tt 


(k—Boltzmann constant, e—electron charge, \—electric conductivity). « was taken 
as 0-5 cal/degC cm/s. As the change of A in metals is inversely proportional to 
the temperature, thermal conductivity, x, of the core does not depend upon 7. 


* If the period of the Earth’s formation including the age of the Earth is equal to about 4-7.10° 
years. 
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The radiative transfer must have a very small effect in the Earth’s core, as the ab- 
sorption of radiative energy in metals is very great:,« reaches 105 cm“. In 
addition the absorption changes very little along the spectrum and does not depend 
consequently upon the temperature and pressure. 

The most uncertain value is the coefficient of absorption e, which is very sen- 
sitive to the increase of T and the radiation frequency. The preliminary value of 
it given by Clark is « = 20-40 cm“ for olivine at room temperature at wavelengths 
of o-8u. As may be judged from the data for semiconductors (Fan 1956) the co- 
efficient of absorption increases with temperature. For T close to 600°K the co- 
efficient of absorption in Si and Ge reaches « = 100 cm~}. 

We took « = 200 cm~! previously (Lubimova 1957) and e = 10 cm~! in the 
present paper. The refractive index was taken to be n? = 3 (about the value for 
olivine or MgO in the visible part of the spectrum). 


6. Conclusion 


Results of temperature calculations and the corresponding thermal conductivity 
distribution is given in Figures 6 and 7. We see that for all the moments of time ther- 
mal conductivity at first decreases with depth, then rapidly increases at great depth. 
The effective conductivity « at the base of the Earth’s mantle is four times larger 
than the ordinary conductivity of rocks if « = 200 cm™! and 15 times if « = 10 
cm-1, Inside the Earth, whence no heat is escaping, the temperature rose rapidly 
with time, while in the first upper hundred kilometres, out of which the heat 
escapes to the surface, a much more gradual temperature increase was taking place. 
Radiative transfer is going on only very deep inside the mantle. The thermal 
conductivity of the upper layers is determined by the molecular component xj, 
owing to which total thermal conductivity « has minimum equal to 0-003 
cal/degC cm/s at a depth of 50-100 km. The value « decreases at this depth by 
15-2 times as compared with the normal thermal conductivity of the surface rocks. 
This decrease of thermal conductivity is followed by a number of peculiarities in 
temperature distribution. On the background of the general decrease of the tem- 
perature gradient with the approximation to the centre a small interval of depths 
up to 20 km (see Figure 6a) may be chosen, where there is an increase of the 
gradient. 

The rapid increase of the temperature gradient over the depth 1000-15000 km 
(Figure 6) is caused by our schematic supposition about the sharp variation of the 
concentration of the radioactive sources. 

The low thermal conductivity of the upper _— prevents considerable heat 
escape from the inside of the Earth. As a result the heating of the interior of the 
Earth is continuing up to the present moment. Thus the conclusion reached earlier 
is confirmed: the region of heat escape to the surface and the internal part which is 
heated may as before be separated. Although the thermal conductivity of the 
internal layers of the Earth is considerably higher than that of the surface layers 
(it equals 0-03 cal/deg C cm s at a depth of 1000 km if « = 10 cm~) no consider- 
able heat escape takes place from the depths exceeding 500 km, as the temperature 
gradient there is much less than the gradient in the upper layers. 

The change of the heat flow with time corresponding to the distribution of 
temperature in Figure 6 is given in Figure 8. 

The straight line signifies the moment of differentiation. The maximum 
value of the heat flow reached 3-4.10~® cal/em?s and took place 2-5.109 years 
ago. The modern value equals 2-3.10~® cal/cm?s, which exceeds twice the 
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Fic. 6.—Temperature distribution in the Earth for different moments 
of time with account of the variations of thermal conductivity with depth. 








™ 
1000+ ° 
9 b 
@ 
a 
s 
3 
a 
£ 
eo 
ia 
500+ 
. iH T ¥ — 
20 40 60 80 100 
Depth: km 


Fic. 6a.—Temperature in the upper layer of the Earth with account of 

variable thermal conductivity «(r, t). a—from a study of thermal history 

of the Earth. b—extrapolation of the observed gradient with account of 
heat generation changes in granite, basalt and ultrabasalt layers. 





ob 
an 
wi 
th 
Ez 





Thermal history of the Earth [129 


observed mean value 1-2-10~ cal/em?s. This value was recently obtained by Clark 
and Niblett (1956) for the region of Alpine railroad tunnels. The heat flow obtained 
will be close to the mean observed for half the heat generation in the crust such as 
that given in Figure 5. This leads to a suggestion that the uppermost layer of the 
Earth cannot consist of a continuous granite massive only; the content of granites 
must be less. It agrees with Birch’s (1955) idea about mosaic structure of the crust. 
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Fic. 7.—Distribution of thermal conductivity in the Earth’s mantle. 


The selection of the correct heat generation in the uppermost layers of the Earth 
is connected with the problem of granite distribution. 

The depth of the region of heat escape is about 500 km. The temperature in- 
side this has decreased about 100°C in the last 10° years. 

We shall compare the temperature at various depths with Uffen’s curve of 
fusion (see Figures 6 and 5). The temperature of the Earth increases rapidly in 
the upper layers after which its rise is lessened. 

Contrariwise, the fusion temperature of the Earth rises slowly in the upper 
layers, in which the pressure is not large and rapidly in the inner layers. This 
leads to a tendency to fusion in the upper part of the Earth’s mantle. This tendency 
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can also be noticed in the curves of Figure 5. The molten stratum formed in limits 
from 100-200 up to 500-700 km. Some suggestions may be made as to consequences 
of the appearance of the stratum. This melting might have been the cause of the 
extensive magmatic activity in the past which, according to certain ideas, led to 
the formation of the modern Earth’s crust. 

Differentiation of matter in the stratum under consideration must have resulted 
in reduction of heat sources in the lower part of the stratum and in their increase 


uw 





Heat flow: 10~° cal/cm? sec. 


9 H 2 = ae? ee 
t: 10° years 
Fic. 8.—Heat flow variation with time after the redistribution of heat 
sources at the time t) = 1-6.10° years. 


at the top of it; it may be connected with the nature of the C-layer and the layer of 
low u, and u, velocities in asthenosphere. Movement of sources towards the surface 
must have led to decrease of temperature gradient at the base of the stratum 
and its increase at the top of this, which would explain the observed increase of 
u, and u, velocities and density p in the C-layer as well as decrease u,, u, in astheno- 
sphere at depth about 150 km. Some accumulation of iron at the lower part of 
the stratum due to the differentiation must have raised the density too. Moreover, 
it may explain the observed rise of electrical conductivity in the C-layer. 

It is feasible that the modern magmatic reservoirs are individually preserved 
oases of melting at the upper boundary of this layer. As to the lower boundary, as 
it is seen from Figures 5 and 6, the zone of the reservoirs is displaced with time 
towards the inside of the Earth. Its limiting depth is 500~700 km for the present. 
This is also the limiting depth for deep focus earthquakes. 

The Earth as a whole could not melt in the process of its evolution. Complete 
melting is only possible for bodies of much lesser mass than the mass of the Earth, 
such as the Moon for example (Levin & Lubimova 1955). However, in the case 
of smaller bodies like the asteroids the melting could take place only in the centre 
(Lubimova & Starkova 1954). If the melting temperature for the metallic state 
of silicates is supposed to be close to that of iron, as estimated by Simon (1953) with 
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Gilvarry’s correction (1957) for the pressure in the Earth’s centre, then, as is seen 
from Figure 6, the outer parts of the core must be in a molten state, while the 
internal core must be solid. This agrees with the hypothesis of Jacobs (1954) 
about the rigidity of the inner core. 

The general conclusion of this paper is as follows. The wide distribution of 
radioactive elements, which is difficult to deny, especially with reference to potas- 
sium, and the small thermal conductivity of the Earth’s upper layers which 
can only decrease with the rise of temperature, must result in heat accumulation 
in the Earth’s interior. But the rise of temperature proceeded so slowly that no 
catastrophe could occur and it only gave rise to slight expansion of the globe. 
The increment of the Earth’s radius is about 7 cm for 1,000 years. 

It is significant that the outer layers were not heated but rather cooled down 
during the last (1-2).109 years. The resulting compression of the surface layers 
covering the expanding bedding may be one of the causes of the formation of elastic 
strains in the upper part of the mantle which in turn result in faults and earth- 
quakes. The possible appearance of a molten stratum in the upper part of the 
Earth’s mantle below the crust gives one possibility of explaining the non- 
uniformity of this part and the origin of magmatic reservoirs. 
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Note added in proof.—An article by Sydney P. Clark on ‘“‘Radiative transfer in 
the Earth’s mantle” appeared in Trans. Amer. Geophys. Union, 38, No. 6, 1957, 
while I was reading the proofs of the present paper. In that article Dr Clark 
comes to the conclusion that ‘‘the presence of an outermost 600 km shell having a 
high conductivity (due to radiative transfer) may increase the heat loss from the 
Earth”. He does not take into account, however, that thermal conductivity of the 
upper first hundred kilometers is determined by its molecular component, which 
decreases as the temperature rises and thereby prevents considerable heat losses 
in spite of high conductivity of the Earth’s interior. 

1958 February. 
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Seismic Observations from the 1956 Atomic 
Explosions in Australia 


B. A. Bolt*, H. A. Doyle}, and D. J. Sutton{ 


(Received 1957 November 19) 


Summary 


A discussion is given of seismic observations from the four 1956 
British atomic explosions, at Maralinga, in South Australia. Both P 
and S phases were observed at field stations set up along the Trans- 
Australian railway for 0-4° < A < 11°, and on a Benioff seismograph 
located near Adelaide at 8°. 

The travel-times indicate a Py velocity of 8-21+0-005 km/s and 
an S, velocity of 4-75+0-01 km/s. There is also evidence for P and 
S phases travelling near the surface with velocities of 6-03 + 0-09 km/s 
and 3-55+0-04 km/s respectively. No onsets corresponding to paths 
through intermediate layers were observed. 

The observations can be explained by a single crustal layer having 
constant seismic velocities. This hypothesis gives a crustal thickness of 
32 km from P and 39 km from S. 

The paper contains a more detailed treatment of the Maralinga 
data previously discussed by Doyle. These data have provided the first 
published studies of Australian crustal structure based on refraction 
methods. 


1. Introduction 


Early in 1956 it was announced that a number of British atomic bombs would 
be exploded at the testing range, near Maralinga, South Australia. Useful seis- 
mological information had previously been obtained from atomic tests in the 
U.S.A., the Pacific and Australia so, at the suggestion of Sir Edward Bullard and 
with the co-operation of the British Atomic Weapons Research Establishment, 
a programme was developed in Australia to obtain further seismic recordings from 
these explosions. 

The programme was in three parts. The first part, undertaken by the Depart- 
ment of Geophysics, Australian National University, consisted of establishing 
field stations along the Trans-Australian Railway, west of Maralinga. Sir Edward 
Bullard originally suggested investigating the 20° discontinuity, but as the chances 
of recording to this distance were small, Professor J. C. Jaeger suggested recording 
along the railway towards Perth (A = 15°). Those involved in establishing and 
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operating these stations were H. A. Doyle, R. Green, H. Angus and E. Irving 
(all of the Australian National University), and B. A. Bolt (of the University of 
Sydney). As the second part, the Bureau of Mineral Resources undertook more 
detailed seismic recording within about 1° of the explosions. The third part 
involved the five permanent Australian seismological observatories, Brisbane, 
Riverview, Melbourne, Adelaide and Perth. Their positions relative to Maralinga 
may be seen on Figure 1. Arrangements were made to inform the Observatory 
directors of the explosion times in advance, through Professor K. E. Bullen. 
A Benioff seismograph was specially installed at Adelaide, by Dr Sutton, to replace 
the Milne-Shaw instrument and Dr J. Webb of the University of Queensland set 
up a field station at Quilpie (Qld.). 

Atomic explosions occurred on September 27, October 4, October 11, and 
October 21, 1956; these will be called explosions A, B, C and D. Seismic waves 
were recorded at the seven field stations shown in Figure 1, and by the Benioff 
instrument at Adelaide, but none were recorded at Brisbane, Quilpie, Riverview, 
Melbourne or Perth. 
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Fic. 1.—Maralinga and the recording stations. 


In this paper, the observations from the field stations and Adelaide are 
discussed; those of the Bureau of Mineral Resources are not included. A prelimi- 
nary study of the field observations has been given by Doyle (1957). Prior to these 
works no estimates of seismic velocities in Australia, comparable with studies in 
Europe, Japan, North America, South Africa and elsewhere, have been published. 
Enough detail is therefore included in the present paper to allow the uncertainties 
to be independently assessed. 


2. The field stations 


Four high-speed recorders and a Willmore seismograph were used. The 
high-speed instruments were built at the Bernard Price Institute, Johannesburg, 
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for seismic work in South Africa, and were kindly made available by the Institute 
for this experiment. They were used with a vertical and a horizontal moving-coil 
geophone (the Forrest instrument had an extra vertical geophone), of natural 
periods about o-6s. These were connected through battery-powered ampli- 
fiers to a galvanometer recorder, with paper speed 10-20mm/s. A third 
galvanometer was connected to a chronometer which marked half-seconds and 
minutes. 

Accurate clock corrections were obtained by recording time signal broadcasts 
on one of the galvanometers before and after the expected explosions. Time 
signals from WWVH (Honolulu) on 15 or 10 Mc/s, and VHP (Belconnen, Mt. 
Stromlo) on 6-4 and 8-4 Mc/s and 44 kc/s were found the most useful. The 
cessation and return of the tone on WWVH could often be read on the record 
within 0-05 s even when the pips were obscured by static. Dr A. R. Hogg, the 
Acting Commonwealth Astronomer, informed us that, because of travel delay and 
errors in adopted longitudes, WWVH signals would be about 0-07s slow 
compared with Belconnen. Since the times of the explosion were known to only 
o-1s this correction was ignored. 

The Willmore instrument was used with a vertical and a horizontal geophone, 
of natural periods 1 and 13s respectively. Minute marks from a chronometer 
were identified by means of hand pulses coinciding with the radio signals. As a 
check, chronometers at all stations were rated by ear each day. This check showed 
that the rate of the Willmore chronometer was not constant during the experiment. 
Corrections from a graph allowed times to be found within o-2s, which is less 
precise by a factor of four than the times from the fast recorders. 

The field stations were established in buildings, kindly made available by the 
Commonwealth Railways, in the settlements along the Trans-Australian line 
The railway telephone system provided convenient communication between 
stations. The Zanthus, Celebration and Southern Cross sites were situated directly 
on the Yilgarn Pre-Cambrian Shield, while the remainder were on the Nullarbor 
Plain. This plain consists of horizontally bedded Miocene limestones which, 
from the evidence of boreholes, are probably no more than a thousand feet 
thick. 

After explosion A, during which five field stations were operating, the number 
of observers was reduced. For subsequent explosions, the remaining observers 


either recorded at stations at which satisfactory results had not been obtained or 
established new stations. 


3- The Mt. Bonython (Adelaide) Station 


A three component, Benioff short-period moving-coil seismograph, with film 
recorder, was used. Initially the seismograph was located on piers in the basement 
of the observatory in the grounds of the University of Adelaide. However because 
of local noise it was moved for explosion D to the concrete floor of a shed at Mount 
Bonython, 11 miles from the city. 

The stratigraphical sequence for Mt. Bonython indicates the presence of up to 
6000 feet of consolidated sedimentary rocks below the site; the height above sea. 
level is 2 200 feet. 

Time pips at half-minute intervals, obtained from the South Australian 
Standard Time Service, were transmitted to the site by land-line. These pips 
differed from Belconnen time by less than 50 ms. 
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4. Observations 


Details of the records obtained are shown in Table 1. The paper speeds are 
given only as orders of magnitude. 


Table 1 
Particulars of Records 
Station Instrument Components Explosion Paper Speed 

mm/s. 
11 Mile Camp B.P.I.(1) Z, H(2) Cc 10 
Cook a Z A 20 
Forrest as Z B 10 
Loongana Willmore Z, H(2) A I 
Zanthus BP 4. Z, H(2) A 20 
Adelaide Benioff Z, 2H D 03 
Celebration BP. I. Z, H(2) C.D 10 
Southern Cross “ Z, H(2) D 20 


(1) High-speed recorder built by Bernard Price Institute. 
(z) The horizontal geophones were orientated roughly towards Maralinga. 


The co-ordinates and origin-times of the four explosions were kindly supplied 
by the Atomic Weapons Research Establishment. 

The nearest recording was obtained at 11 Mile Camp using Explosion C. 
The record is very clear (Figure 3) but in this test the atomic device was dropped 
from an aircraft, causing higher uncertainty than usual in the location and time of 
the explosion. The supplied angular distance and explosion time have accuracies 
of +0°-o1 and +18, respectively. The combined error for one observation is, 
however, unlikely to exceed 1-5s which is the usual value of the standard error of 
reliable P observations from earthquake work (Jeffreys 1947), and the observations 
were useful for comparison with the calculated travel-times. 

The co-ordinates and times of the remaining explosions have accuracies within 
+0-1 minutes of arc and +0-1s respectively. The positions of the field stations 
were obtained from railway survey data and Department of Lands maps, and their 
accuracy is +0°-005. Angular distances A were calculated from the geocentric 
co-ordinates with six-figure tables using 2(1—cos A). Bullen (1934) has shown that 
the resultant maximum error is 0°-034 for 1° < A < 20°. The full ellipticity and 
height corrections are less than 0-05 s for all observations and were ignored. 

Each author independently examined the records, reading only the prominent 
phases, and calculating the time corrections. On comparison, the three sets of 
travel-times were in striking agreement, only a few readings differing by more than 
o-2s. Those readings supported by only one of the authors were omitted; the 
remainder are shown in Table 2. 
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Table 2 
Arrival Times of Phases 


Station A Travel-time Phase Weight Residual 
deg. s s 
11 Mile Camp 0°47 8-2 a (3)* —0°45 
14°9 S, (3) +020 
I51°3 Sound 3 
wave 
Cook 1°32 24°0 P, 3 —o'26 
25°7 P, 3 +0°59 [ +0°03] 
41°9 Si 3 +068 
44°9 Sn 2 —0'24 
60-0 2 
Forrest 3°19 50°4 | 3 —o'03 [—0-26] 
58-0 P, 3 —o-76 
62°1 3 
89°5 Sn 3 +0°54 
97°6 S, 3 —2°24 
Loongana 4°13 63:2 Pe 3 +0-04 [—0-12] 
72°0 2 
82-0 P32 I +5°77 
II1‘2 Sn 3 +004 
131°5 S; 2 +1°97 
Zanthus 7:08 103'I m 3 —o-'00 [+0-41] 
110°8 2 
132°1 P, 2 +1°53 
179°9 Se 3 —0'29 
222°4 S, 2 +0°52 
Mt. Bonython 7°85 116-2 P. 3 +2°67 [+3:°14] 
203°5 S. 3 +5:°22 
213°0 2 
248°5 S, * +2°29 
Celebration 8-72 
Explosion D 125°3 P, 3 —o-o1 [+0°73] 
219°4 Sn 3 +0°45 
272°0 S; 3 —1-°20 
Explosion C + 126°5 igs (1) +119 [+1-93] 
217°8 ae 2 —0'75 
Southern Cross 10°69 150°2 Fs 3 —1°88 [—o0-76] 
264°4 Sn 3 —0°43 


*() denotes observations, which although weighted, have not been included in the formal 
calculations. 


A simple system of weighting has been adopted based on the degree of agree- 
ment in the selection of phase. Weight 3 is given to all onsets agreed on within 
0-2s by the three authors; weight 2 is allotted when only two authors agree 
within 0-2; weight 1 indicates a spread of up to 1s in the independently adopted 
arrival time of the phase. 
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The records of the fast recorders were the simplest to interpret, showing little 
complexity after the first arrivals of the main phases. This is probably due to the 
shorter period of the instruments. The vertical components had the larger ampli- 
tudes and the sharper onsets but the horizontal component provided a useful 
check. The Willmore records had many secondary onsets with often no agreement 
between vertical and horizontal components. This has also been commented on 
by Hodgson (1953). The Benioff instrument gave clear records showing small 
amplitudes, the largest being on the vertical component. (Although the first P 
onset was small the reading is not in doubt.) 
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Fic. 2.—Travel-times for all stations. 


The travel-times are shown plotted against distance in Figure 2. This strongly 
indicates that there are four phases with characteristic velocities over the profile. 
These phases have been called Pn, Pi, Sn, Si in the order which corresponds to 
the order of arrival beyond Cook. Examples of the phases are reproduced in 
Figure 3. 

In a similar way to the observations from the Heligoland explosion, reported 
by Willmore (1949), the first arrivals are uniformly well recorded over the whole 
range. In this respect the present observations are also in agreement with Tatel 
and Tuve’s findings from blasts in Tennessee and in disagreement with Gutenberg’s 
suggestion, based on Californian earthquake studies, of a shadow zone for P», 
for 7° < A < 13° (cf. Byerly 1956). 

The Heligoland S observations are, however, dissimilar. Willmore states that 
transverse waves were poorly recorded by the short period vertical instruments; 
in the Maralinga experiment the S phases were generally comparable to the P 
in amplitude but had slightly longer periods (0-1—0-2s). Gutenberg (1946, 
1952) has drawn attention to S phases recorded from American atomic tests. 
From the 1945 New Mexico test they were recorded with indistinct beginnings 
and small amplitudes at observatories having 3°-9 < A < 9°-7; from the Nevada 
tests they were recorded up to A = 3°-6. In the latter case Gutenberg remarks 
that it is possible that some of the S phases are surface waves. 
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5. Travel-times 


Least squares solutions for the travel-times of P,, P;, Sn, Sj were made using 
linear equations and allotting weights as indicated in Table 2. Equations for the 
Mt. Bonython and 11 Mile Camp observations were not included since Mt, 
Bonython is in roughly the opposite azimuth to the field stations and because of the 
relatively large uncertainty in the angular distance between explosion C and 
11 Mile Camp. 


For Pp, using only the first arrivals at Forrest, Loongana, Zanthus and Cele- 
bration, the result is: 


t = (7°24 + 0°05)s+(13°54 + 0-01) s/deg. 


The P, surface velocity is (8-21+0-005) km/s and the standard error of one 
observation o = 0-048. This gives residuals of +0-59 and —1-88s for Cook and 
Southern Cross respectively. The Cook residual is quite satisfactory considering 
that the phase is a second arrival, but the Southern Cross residual is rather large 
and suggests the inclusion of a cubic term in the travel-time equation. An alter- 
native solution, obtained by including equations for Cook and Southern Cross 
gives: 
t = (8-10+ 0°49) s+(13°36 + 0°07) s/deg. 


The corresponding velocity is (8-32 + 0-04) km/s o and increases to 0-58s which 
is somewhat high for clear P, arrivals measured independently by the three authors 
to within 0-2s. The residuals from this alternative solution are included in Table 2 
in brackets. 

For distances up to 14° the J.B. times include a cubic term, approximately 
—o-0025(A—1)8, which amounts to —1-14s at 8°-7 (Celebration) and —2-3s 
at 10°-7 (Southern Cross). These values may be compared with the residuals 
from the first solution, —o-o1 and — 1-88s respectively. Jeffreys (1952) has found 
that the times in Japanese earthquakes agree reasonably well with the J.B. cubic 
term, while times in Europe, and less certainly California, are consistent with 
linear formulae. 

Observations marked P, in Table 2, from Cook, Forrest, Loongana and Zanthus 
give t = 0-728+18-43A which makes the Loongana residual +5-1s. The 
complexity of the Willmore record after the first arrival makes the identification 
somewhat uncertain so the Loongana equation was omitted, giving t = o-gos 
+ 18-71A with o = 0-88s. The constant term therefore is not significant and we 
take it to be zero. This is equivalent to assuming a single-layered crust having 
constant P velocity. The solution is 


t = (18-45 + 0-28) s/deg. 


The P, residuals in Table 2 are obtained using this formula; o = 1-49s. The 
surface velocity is (6-03 +.0-09) km/s. Support for this solution comes from the 
first arrival at 11 Mile Camp, which has a residual of —0-45s. We note that on this 
solution P, would not appear before 75 km, so the second arrival at 11 Mile Camp 
is not Pp. 

There is no difficulty in combining the S, observations from Cook, Forrest, 
Loongana, Zanthus, Celebration and Southern Cross. The solution is 


t = (14°33 + 0°46)s + (23-43 + 0°06) s/deg, 


corresponding to a velocity of (4-75 +0-01) km/s over the surface. The standard 
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error for one S, time is 0-48s. As with P,, the residual at Southern Cross is 
negative (—0-4s). The J.B. formula for S, at the distances involved here (Jeffreys, 
1939) contains —0-00314(A—1)° which contributes — 31s at 11° and is therefore 
inconsistent with the Maralinga data. 

The five S; observations at Cook, Forrest, Loongana, Zanthus and Celebration 
provide a check on our hypothesis of a uniform single-layered crust. We obtain 


t = (1-764 1°94)8+(31°00 + 0°36) s/deg. 


with o = 1-97s. This standard deviation is of the order of the standard error of a 
good S, observation from earthquakes (Jeffreys 1939), but unexpectedly high for 
astudy with known distances and origin times. Until further information becomes 
available it seems best to suppose the observations are the same phase. The con- 
stant term may then be dropped to give 


t = (31°35 +.0°40)s/deg, 


corresponding to a surface velocity (3-55 +0-04) km/s. The residuals against this 
solution are given in Table 2; o = 1-80s. The reduction in o is due to reduction 
in the large unexplained residual at Forrest. The residual of +0-2s for the 
second arrival at 11 Mile Camp agrees with the solution and we conclude that the 
S; data is consistent with a single-layered crust. 

The Mt. Bonython observations now give the following residuals, P, + 2-78, 
Sn +5-28, and S; +2-3s. Delays of this magnitude are rather surprising. Explana- 
tion of the P, and S, residuals in terms of large down-faulted blocks in the Adelaide 
“shatter zone” leaves the Sj residual unexplained. It seems best to wait for further 
data from Mt. Bonython for a final answer. 


6. Discussion 


The velocities and adopted crustal structure give the following crustal thick- 
nesses : 
From P waves 32-2 km. 
From S waves 38-8 km. 


The various standard errors imply uncertainties of order 3 km for both estimates. 
Hodgson (1953) and others have similarly obtained a higher value for the crustal 
thickness using S. 

The phases P; and Sj are of particular interest. Identification of P; is some- 
what in doubt at stations beyond Forrest, while S; is clear at A = 8°-72 (Celebra- 
tion) and its presence on the Southern Cross record may have been obscured by the 
high noise level at the time. In their study of the South Australian earthquake, 
1939 March 26, situated about 10° east of Maralinga, Bullen & Bolt (1956) identified 
two prominent phases on the Adelaide record. These arrived 12s and 48s after 
P,, and were tentatively called P* and S*. The P;— Pn, Si— Sp intervals from the 
Maralinga travel-time curves at the same distance, A = 3°-0, are 8s and 47s 
respectively. The agreement between S* and 5; is an interesting feature since the 
S* was thought to travel in an intermediate layer. 

In this paper the P; and S; phases are taken to be waves travelling close to the 
surface. There is some indication that they correspond to the longitudinal and 
transverse waves II, and L,. The L, phase has been found to be efficiently pro- 
pagated in all continental crusts with a mean velocity close to that found in North 











144 B. A. Bolt, H. A. Doyle and D. J. Dutton 


America (Press & Ewing 1952), namely (3-51+0-07) km/s. For Australia, Bolt 
(1957) has found a mean L, velocity of (3-50 + 0-07) km/s. The II, phase has been 
observed by Press and Gutenberg (1956) on seismograms of a North American 
earthquake; it is characterised by small amplitudes and erratic occurrence. They 
obtain a mean velocity of (6-09 + 0-06) km/s from 11 observations with 5° < A 
< 3°. 

There are a few onsets listed in Table 2 which are not explained. The second 
P phase at Loongana may be PP, while the second S phase at Mt. Bonython may be 
SS. 

A large number of studies of travel-times from blasts in other countries are 
available for comparison with the Maralinga results. Of special interest is their 
agreement, set out below, with the results of Hodgson (1953) in the Canadian 
Shield, and with the results of Willmore, Hales & Gane (1952) in the Western 
Transvaal. 


Canadian Shield Western Transvaal Western Australia 
P,, (km/s) 8-18+0-01 8-27 +0-04 8-21 +0-005 
P, (km/s) 6:25 +0°02 6-09 +0°03 6:03 +0°09 
S, (km/s) 4°85 +0°10 4°83 10-02 4°75 to-o1 
S, (km/s) 3°54 10°02 3°68 +0-02 3°55 40°04 
Depth to Mohorovici¢ 
discontinuity (km) 36 36 36 


It would seem that the old continental shields are seismically very similar. 


7. Air waves 


The accuracy of sound velocity measurements over a long distance in air is 
seriously affected by variations in temperature, wind velocity and humidity. 
From experiments over distances up to 14 km, E. Esclangon (1919) gives 339-8 + 0-1 
m/s for the velocity of sound in dry still air at 15°C. Lately, some observations of 
air waves recorded from atomic tests have been published. Gutenberg (1946) lists 
the velocities of direct air waves recorded at five seismological stations, having 
4° < A < 10°, from the New Mexico test. The mean velocity is 350+4 m/s. 
He gives no similar figures for the Nevada test (Gutenberg, 1952) but states that 
many records indicated a phase with a velocity of slightly over 300 m/s. 

Doyle (1954) recorded air waves from the 1953 atomic explosions in South 
Australia. At 18 km a velocity of 340+4 m/s was obtained. No temperature 
corrections are available. 

Air waves from the 1954 hydrogen bomb explosions, recorded in Shionomisaki, 
Japan, at about 36°, have been studied by Yamamoto (1955). In a rediscussion 
Burke-Gaffney & Bullen (1957) give a mean velocity of 321 m/s with an uncertainty 
of 5 m/s. 

In the Maralinga tests a direct air wave was recorded at 11 Mile Camp from 
explosion C. The recording is extremely complex but the first onset can be read 
within o-1s. The travel-time is 151-3s and the corresponding velocity 344 m/s. 
The uncertainties for explosion C (discussed in Section 4), give an uncertainty 
of about 6 m/s. Since the various meteorological conditions are not known, the 
result is in reasonable agreement with previous observations. 
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The Motion of the Surface of a Uniform Elastic Half-Space 
Produced by a Buried Torque-Pulse* 


C. L. Pekeris and I. M. Longman 


(Received 1958 January 23) 


Summary 


In this investigation we determine the motion of the surface of a 
uniform elastic half-space produced by the application of a torque-pulse 
at a point below the surface. The axis of the torque is taken to be parallel 
to the surface and its time variation is assumed to be represented by the 
Heaviside unit function H(t). Our results are compared with those of 
Pinney (1954), who treated the same problem by a different method. 
The principal feature of physical interest which we found is that, for 
ranges where the direct S wave is preceded by a diffracted SP wave, 
the displacement at the surface starts with an infinite amplitude at the 
time of arrival of the SP wave, and that this is followed by an even 
stronger infinity at the time of arrival of the S wave. Also, for small 
ranges, for which there is no SP wave, the displacement starts with a 
sharp pulse in the form of a Dirac delta function. None of these 
features was brought out by Pinney’s curves. The results are shown in 
Figures 2, 3 and 4. 


1. Introduction 


In recent investigations on the problem of the propagation of a pulse of the 
form of a Heaviside unit function H(t) in doubly refracting media, the following 
result of physical interest was obtained. When the slow-speed wave is incident 
on an interface at an angle greater than the angle for total reflection, it is partially 
transformed by diffraction into a high-speed wave, which is propagated along the 
boundary before re-entering the medium. This high-speed diffracted wave arrives 
ahead of the direct slow-speed wave, even though its trajectory is longer. The 
exact solution of the wave equation for this pulse-problem shows that the amplitude 
becomes logarithmically infinite at the time of arrival of the reflected slow-speed 
wave. On the other hand, no such infinity exists for shorter ranges from the source, 
which are inside the critical distance for total reflection. This was found to hold 
for electromagnetic (Pekeris & Alterman 1957), acoustic (Arons & Yennie 1950, 
Pekeris & Longman 1958, Jeffreys & Lapwood 1957) and seismic (Friedlander 
1948, Pekeris & Lifson 1957) pulses. In the acoustic case some experimental 
evidence in support of this theoretical result was produced by Arons & Yennie 
(1950). It was also suggested (Pekeris & Longman 1958) that this effect may have 
a bearing on the observations of Tatel & Tuve (1955) on the behaviour of the 


* This work was supported by the Office of Naval Research under Contract Nonr—2189(00). 
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reflection from the Mohorovici¢ discontinuity in the Earth’s crust as one passes 
the critical range. ; 

In this investigation we wish to determine the motion of the surface of a uniform 
elastic half-space produced by the application of a torque-pulse at a point situated 
at a depth H below the surface. The torque is taken to be oriented parallel to the 
surface, with its axis pointing in the positive x-direction, as is shown in Figure 1. 
The time variation of the torque is assumed to be represented by the Heaviside 
unit function H(t). The source emits shear waves (S) only. On reaching the 
surface, the incident S wave is reflected partially as an S wave and partially as a 
P wave, if the angle of incidence is less than the angle for total reflection. For 
ranges beyond the angle of total reflection, we have from geometrical optics that 
all of the incident energy is reflected as an S wave. However, because of the 
sphericity of the incident wave front, a diffracted P wave is propagated along the 
surface. We shall designate this diffracted wave as the SP wave. Under these 
circumstances a point on the surface first receives the SP wave and then the direct 
S wave, although the path of the former is longer. 
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Fic. 1.—A torque-source, whose axis is oriented along the x-axis is applied 
at the point O below the surface. @ is the angle in the x-y plane; ¢ is the 
angle in the y-z plane. 


In a previous investigation (Pekeris & Lifson 1957) of the case of a buried 
seismic point source of the form of a downward force varying like H(t), it was found 
that at the time of arrival of the S wave the amplitude becomes logarithmically 
infinite. For the problem at hand of the buried torque-source, one would expect 
an even stronger infinity for the displacement at the surface. Now the buried 
torque-source problem was treated previously by Pinney (1954), and his computed 
curves for the displacement of the surface do not show any infinities. We have 
therefore undertaken to resolve this problem by a different method, and to compare 
the results with Pinney’s solution. It turns out that both solutions agree and give 
strong infinities at the epoch referred to, as well as at other times. The reason 
why Pinney missed these singularities, which are the features of the solution 
of principal interest to the seismologist, is that his method of numerical evaluation 
of the integrals appearing in his solution becomes inaccurate just where the 
singularities occur. 


2. Formal solution 


We adopt the system of coordinates shown in Figure 1, with the x-axis 
K 
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pointing in the direction of the axis of the applied torque. Let the displacement 
vector u be given by 


u = grad¢+curlT, (1) 
then the boundary conditions can be satisfied by taking 

Tz = Q, Ty = 0, Tz = OP/dx, d = Op/dy. (2) 

The components of u, and the dilatation A are then given by 
u = $2y+ Pry, (3) 
v = Pyyt+ Q2—Prz, (4) 
w = pyz— Qy, (5) 
A = divu = Vy, (6) 


where the subscripts denote partial differentiation. On substituting into the 
equations of motion, we find that they are satisfied if 


V2Q—k2Q = 0, (7) 
V2P—RP=0, V2sb—hb = 0, (8) 

where 
He = ple, 1H = peg, c= lp, cg? = Q+2u)p= 3% (9) 


Here c denotes the velocity of shear (S) waves, cp the velocity of compressional 
(P) waves, and p denotes the operator 0/d0t. A and yw denote the Lamé elastic 
constants of the medium, which are assumed to be equal. 

Consider now the solution of the wave equation (7): 


Q = R'f[t-(Ri/c)], R= (x2+y24+ 22), (10) 


where f is an arbitrary function. 
Taking a cylindrical system of coordinates (r, ¢, x), as shown in Figure 1, 
with r = (y? +22)#, we find from (1) and (2) that O represents a torsional wave 


with only one component of displacement u, along circles in planes normal to 
the x-axis: 


_- - «- =. oo ce 





0Q Ou, a2Q pre 0 (1 OQ (1) 
- = = — — }. II 
( R OR 
Let T’ denote the torque which the matter to the left of a plane normal to the 
x-axis exerts on the matter to the right of this plane. We have 


oo 


2f1R=° 
T’ — al lea = anus |(RE—20)( F an 7) Rut 


= amus{-A-) + =flt—(aio)]}. (12) 
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In the limit of vanishing x we get 


T’ > 4muft). (13) 
The same torque 7” is exerted at O by the matter on the right of the yz plane on 
the matter to the left, so that the total torque T exerted by the source represented 
by Q in (10), is 

T = 8nufit). (14) 
The operational form of the solution of the wave equation (7) for a torque-source 
of magnitude JT and varying with time like H(t) is therefore given by 





Qo 


a “= ( t 5) 


One of the boundary conditions is the vanishing of all the components of 
stress on the surface, which is a horizontal plane: 


Teg = Tyg Ff Teg = Oo, z= —H. (16) 


Accordingly we use the representation 





Q = T e-kR _ fo f Eeaney (Er) dé (17) 
SauR Sm JAB , 

with 
r= (x®+y?)t, RB = (€2+h?)t, kaw = [£2 + (h?/3)]}. (18) 


Following the method of Lamb (1904), we represent the solution of the wave 


equations (7) and (8) by é-integrals like (17), and, denoting the integrands by 
bars, take 


yp = Ae-2Jo(ér), P= Be-*éz] (Er), (19) 
Q = (T/8mp)(E/kB)e***\Jo( Er) + Ce-#P2Jo(Er). (20) 
The three boundary conditions (16) yield 
A = (T)/8mpM)4E(2€? + h2)e-HerOH, (21) 
B = —(T/82pM)8Eék2aBe-2kFH, (22) 
C = —(7/8mpM)(E/kB)[(2k? + k?)? + 47h? o Be 2kAH, (23) 
with 
M = (262+ #)2— 42h. (24) 


Substituting in (5), and putting = kx, we get for the vertical displacement wo 
at the surface 


es TR sin 6 [ x?[x? + (4) ]#exp[ — RH (x? + 1)*]Ji(Rrx) dx 
{(2x2+ 1)?— 4x2%(a2+ 1)ifx2+ (4) 





(25) 


27h 
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where @ is measured from the direction of the x-axis and sin 6 =y/r (see Figure 1), 
Similar expressions can be derived for the other components of displacement, 
In the sequel we confine our attention to wo only. 





3- Operational interpretation 


The operational interpretation of (25) can be readily carried out by the method 
developed by the senior writer (Pekeris 1956; see also Cagniard 1939, Pinney 
1954). It can be expresed in the form 


wo = (T/72uRr) sin 0(0S/0r), (26) 
where 


r=(ct/R), 2 = [2+ B= @2+x), fle) = 


a 
[(2x2 + 1)2— 4220p] 
and S is defined below. The operational interpretation of (25) differs in the two 
cases when the horizontal range 7 on the surface is less than or greater than the 
critical distance H/1/2 for total reflection. In the latter case the motion becomes 
more complex due to the appearance of the diffracted SP wave. 
We have for r < H/,/2: 





(27) 


S = 0, tT <I, (28) 


itr 
a 


S = S = Re | B(r—hB)f(x)dp, = + >t. (29) 


0 


Here Re denotes the real part, and 
B = hr+iposing, — po = [(x—Fr2— 1), (30) 
h=(H/R) x= (8-1), « = [P?—(8)]. (31) 


When r > H/1/2, the first arrival is the diffracted SP wave, at a time 
corresponding to 


T* = (§)'h+(3)'(1 — A), (32) 
while the direct S wave arrives later at r= 1. We have then 
S =o, 7 < 7%, (33) 
tr 
S= —S,= —- im| B(r —hB)f(x)[(«? + sin? 6)-*]cos ddd, 7* < 7 < 1. 
. (34) 


Here Im denotes the imaginary part, and 


5 = [(§)'—Ar], = m = (8 — 4? mm?) em = [(1—F\(1—7*)]#, (35) 


B = hrim(x?+sin?¢)!, x = i(1—f2)t, « = i[(#)—f?]. (36) 
After the arrival of the direct S wave, we get 


S = §,, t > (§)t/h, (38) 
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S3 = 61 | B(r — hB)f(x)[(u0? + 8? sin? ys) cos p] dy, (39) 
0 
B = hr+65 siny, x = 1(1—f?)}, a = 1[(%)—f?]*. (49) 
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Fic. 2.—Vertical displacement at the surface wp» as a function of time. 
W = (167°H?/Tsin 0)wo, 7 = (ct/R), R = (x? +y?+2%)8. 


4. Discussion of results 


We have evaluated the quantity W, which is related to the vertical component 
of displacement at the surface wy by 


1672 pH? 16H2 aS 


W = ——~up = ——_ —. 
T sin8@ . Rr Oar (41) 
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W has also been computed by Pinney (1954) by his method*, thus allowing direct 
comparison of results. The integrals S}, S2 and S; in (29), (34) and (39) were 
computed on WEIZAC using a 48 point Gaussian quadrature formula (Davis 
& Rabinowitz 1956). In order to determine S at a given value of the reduced 
time r = (ct/R), one or two (as in Equation 37) of these integrals had to be evaluated, 
aS/dr was then determined from [.S(7 +e) —S(7 —«)]/2e, using values of 0-001 
and o-ooo1 for e, and the extent to which these results agreed was taken to indicate 
their accuracy. 
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Fic. 3.—Vertical displacement at the surface w» as a function of time. W = (167*H?/T sin 09)», 
+ =(ct/R), R = (x?+y?+27)!, 
SP denotes the diffracted wave; Sthe shear wave. A:r =10H; B:r =16H. 


Figure 2A shows W as a function of 7 for the range r/H = 0-25 for which 
the diffracted SP wave does not yet exist. The quantity (16H?/Rr)S(r), whose 
7-derivative is equal to W, is shown in Figure 2B. It is seen that the latter curve 
starts with a finite jump at tr = 1, which contributes a 5(t) (Sherwood 1958) function 
to W at this epoch. This 8(t) pulse is not shown in Pinney’s curve. Otherwise 
the results agree. 

Figure 3 shows W as a function of 7 for the ranges r = 10H and r = 16H, 
for which the diffracted SP wave exists. The SP wave starts with an infinite 
amplitude. The singularity is of the form —(7—7*)-+. After this sharp arrival, 
the displacement remains small until the arrival of the S wave at tr = 1, when it 
becomes infinite like (r —1)~1. Neither of these singularities is shown in Pinney’s 
curves. 

Figure 4 shows an enlargement of the curves in Figure 3 near the time of 
arrival of the S-wave at tr = 1. The triangles and the dashed portions of the curves 
represent Pinney’s results. 


* Pinney’s 7 is equal to our 7 multiplied by R/H. 
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Fic. 4.—Enlargement of the curves in Figure 3 near r = 1. The triangles and the dashed 
portions of the curves represent Pinney’s results. A:r = 10H; B:r = 16H. 


We have also recomputed W using Pinney’s formulae*, and the results obtained 
are in complete agreement with the values derived by our method, except for the 
delta function referred to above. 
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Summary 


A new set of earthquakes is studied to check the times of P in 
European earthquakes up to 30°, and especially the behaviour near 20°. 
The corrections to the author’s table of 1939 agree closely with those 
found from an earlier set in 1935. Lehmann’s observations of Pd are 
consistent with her interpretation if Pr overtakes Pd at about 14°, but in 
any case dt/dA for Pd is remarkably constant. It seems probable that 
there is triplication and that the observations refer to a neighbourhood 
of the more distant cusp on the time curve. 

No useful new data were available for the Japanese region. 

Several earthquakes in North and Central America showed that 
from 8° to 30° the 1939 table is correct to a fraction of a second apart 
from a possible additive constant. Additional information on the times 
at distances under 8° was given by a paper of Lehmann. The data in 
this indicate that the European velocity is nearly correct for North 
America, but that, relative to 5°, the times from 8° to 30° for American 
earthquakes are about 2s shorter than for European ones. 


1. European earthquakes 


In two former papers (Jeffreys 1952, 1953) I examined new sets of earthquakes 
with “European”, Japanese and American epicentres. Corrections to the 
Jeffreys—Bullen P table (Jeffreys 1939) were found from European and Japanese 
epicentres; the results from American ones were unsatisfactory, partly through 
internal discrepancies and partly on account of the difficulty of separating 
corrections to the table from errors of the epicentres. Substantial corrections, 
mainly owing to a negative correction to dt/dA up to 8°, were found for Europe; 
but those in the specially interesting region about 20° were still doubtful. 

In this paper I study a new selection of European earthquakes from the I.S.S. 
These were chosen because all had observations about 20°. As before, the observa- 
tions were classified by quadrants before analysis, to reduce the effects of errors 
in the epicentres. The useful quadrants are indicated in the third column. 


Precautions as in my papers of 1952 and 1953 were taken to avoid errors due 
to focal depth and multiplicity. 


The observations were grouped by ranges: 2°-8°, 8°-12°, 12°-15°, 15°-17°°5, 
18°, 19°, 20°, 20°-5-25°, 25°-30°. Corrections were expected to vary little in each 


* The first paper on this subject appeared in Mon. Not. R. Astr. Soc. Geophys. Suppl. 6, 348, 
1952. 
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Table 1 
List of earthquakes used 
Date Useful Epicentre Date Useful Epicentre 
Quadrants Quadrants 
1946 Feb. 12 N, E Spain 1948 Feb. 9 N,E,W_ Levant 
Feb. 21 N, W Levant Mar. 6 N Levant 
Apr. 5 N,E,W _ Levant June 30 N,E, W italy 
May 11 §S N. Atlantic July 8 S N. Atlantic 
May 30 N,E,S Switzerland July 24 N, W Levant 
July 16 N, E,W _ Levant Aug. 27. N, W Jugoslavia 
Oct. 13 N, W Levant 
1949 Jan. 14 N, W Turkey 
1947 May 11 N, W Italy Feb. 5 N, W Turkey 
June 4 N, W Levant Mar. 2 S N. Atlantic 
Aug. 30 N, E,W _ Levant May 13 W Turkey 
Oct... ¢ WN Levant July 8 S N. Atlantic 
Oct. 6 N Greece Aug. 17 N, E,W Caucasus 


range except the first, and the I.S.S. residuals (against the 1939 table) were used 
directly. For the first range the observations were grouped at intervals of 1° 
and compared with my previous table (1953, p. 563). This introduced a dis- 
continuity into the standard of comparison, but this could be treated later. The 
data gave standard error 1-84s on 415 degrees of freedom. Comparisons between 
distances by the method used previously gave the corrections in Table 2 (in order 
of estimation), the range 15°-17°-5 being taken as standard. 





Table 2 
Range Correction Weight x? d.f. Entries 
deg s rejected 
15-17°5 ° 
20°5-25 —o'!l 44 30°9 28 ° 
25-30 +0°2 51 34°2 31 ° 
12-15 —0'5 53 34°6 29 ° 
8-12 +0°5 52 28-6 30 5 
2-8 +2°5 42 30°8 25 I 
18 +or1 33 23°2 20 I 
19 —o-6 42 16°1 23 ° 
20 —0'4 30 13°9 21 3 
Totals 212°3 207 


The total y? is normal, but has been obtained as a result of rejecting some entries. 
(An entry is a mean residual over a range of distance for one earthquake.) The 
explanation is partly that the law of error is not normal, and in estimating the 
means and standard error this was dealt with by rejecting residuals over 5s. If 
the law was normal with o = 2s this would make little difference, and previous 
trials have shown that with the actual law it will not habitually underestimate o 
by more than ro per cent or so. Several of the rejected entries depended on a 
single discrepant observation, which would have been rejected at an earlier stage 
if there had been a standard of comparison. What may be more serious is that the 
method compares the mean residual at each stage with the whole of those at earlier 
stages, as if these were exactly known. For instance, if the difference between 
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15°-17°5 and 20°-25° was altered by its standard error, this would alter the 
value for 25°-30° obtained by comparison with the two together. For this reason 
the weights at this and later stages are too high and might be reduced by 30 per 
cent or so. The correct treatment would of course be to set up a least squares 
solution for eight differences of corrections, but this has not been thought worth 
while. ‘The values for 18°, 19° and 20° are all in comparison with the same standard 
and should be comparable between themselves. 


At 20° one rejected entry depended on a single observation. The others rejected 
were as follows: 


Obs. at 20° Weight Final mean x’ Crude 
residual residuals 
Ss Ss 
1946 Oct. 13 N 5 3°4 =—30 9 —3-4-4-2-2 
1948 July 24 W 4 2°7 +4°I 14 o+1+2+41 


The readings are very consistent in both cases, but the means flatly contradict 
each other. In my earlier treatment (1952, p. 354) a mean residual of —1-4s, 
weight 16, was found at this distance, and I remarked that omission of a single 
series (1938 April 19, W) would have changed this to —o-1s. This was what led 
to the doubt at this distance. There was then one discrepant entry; here there 
are two in opposite directions. It is perhaps not out of the question that with 


Table 3 
Mean residuals 
Mean Mean 
Earthquake Quadrant Residual Weight Earthquake Quadrant Residual Weight 
s 8 

1946 Feb. 12 N —o'9 15 1948 Feb. 9 N +0°7 20 
E +09 2 E —0'4 2 
Feb. 21 N —0'5 12 W —I'l 14 
WwW +06 13 Mar. 6 N —o'6 12 
Mar.17 N +11 4 June 30 N —1'2 25 
Apr. 5 N +o'2 16 E +06 6 
E 4+1°7 6 WwW —1'o 8 
WwW +2°1 5 July 8 iS) —0'7 13 
May 11 Ss —O'5 9 July 24 N —o'8 25 
May 30 N —3°0 12 WwW —2°7 8 
E —o'!l 7 Aug. 27 N —0'5 13 
Ss +1°4 6 WwW +0°'9 9 

July 16 N +or'1 17 
E —16 3 1949 Jan. 14 N —o'8 8 
WwW +2°6 8 W —o'9 14 
Oct. 13 N +04 II Feb. 5 N +1°2 8 
Ww —4°7 7 Ww —2°4 7 
1947 May 11 N +o'8 20 Mar. 2 Ss —1'2 17 
WwW —0O'5 15 May 13 W —1°7 20 
June 4 N —2°0 17 June 17 Ss +1°1 17 
WwW —o'6 16 July 8 Ss +1°2 IS 
Aug.30 N —0o'7 25 July 10 WwW +o'8 6 
E +3°1 4 Aug.17. N —1I4 8 
WwW +3°8 6 E —1'2 II 
Oct. 3 N +0°4 10 WwW —o'2 21 

Oct. 6 N —1°4 22 
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a departure from the normal law there might have been a few cases in over 200 
where observations departing seriously from the true value should agree among 
themselves. Whether this is the explanation of the anomaly or not, we have 18 


consistent entries, and can see that these two are in some way exceptional. 


The mean residuals after allowance for the corrections are as follows. They 
may be useful if the epicentres are recalculated. 
To make the results comparable with my 1953 ones, which are adapted to a 
surface focus, 2-5s must be subtracted. 





Table 4 
Comparison of present and 1953 results 
A Present A 1953 
deg. Correction Weight deg. Correction Weight x= df. 
s s 

2-8 oho) 2-8 oho) 

8-12 —2°0 52 8-12°5 —3°7 9 12 20 
12-15 a ng 53 12°5-17°5 —3°6 (20) 25 
I5-17°5 —2°5 (42) 15-19 —2°8 52 49 4! 

18 —2°4 33 19-20 —3°9 30 28 24 

19 —3°1 42 20-25 —2°2 64 43 25 

20 —2'9 30 25-30 —2°9 31 20 13 
20°5-25 —2°6 44 
25-30 —2°3 51 


The bracketed weights are for the comparison of the standard distance with 
2°-8°. 


Table 5 
Smoothed times 
Present 1953 
A Time FY 3? Time r) 5? 
min s s min s s 

14 3 19'0 3 19°0 
13°3 13°7 

15 32°3 —0°3 32°7 —o'9 
13°0 12°8 

16 45°3 —0'4 45°5 —O'4 
12°6 12°4 

17 57°9 —0°3 57°9 —0°3 
123 I2‘1 

18 4 10°2 —0%3 4 10° —0%3 
12°0 11°8 

19 22°2 —o'6 21°8 —0'4 
11°4 11°4 

20 33°6 —0'5 33°2 —0'5 
10°9 10°9 

21 44°5 —o6 44/1 “4 
10°3 10°6 

22 548 —o'2 54°7 —0'4 
10°! 10°2 

23 5 4°9 —o'2 5 4°9 —0'2 

9°9 10°0 
24 14°8 —o'l 14°9 —0'4 
98 9°6 
25 24°6 24°5 











158 Harold Jeffreys 


The contributions to x2 show no inconsistency, though the data were 
combined in different orders and somewhat different ranges of distance. The 
principal difference is about 20°, —o-gs, weight 15, corresponding to a standard 
error of about o-5s and it is not remarkable for an extreme value. The present 
treatment goes a little more into detail than the previous one (but some of the 
information was recovered in the previous one by a reclassification of the residuals 
at shorter intervals). The previous corrections have been smoothed (Table V, 
p- 563) and the only change worth making is to increase the time at 20° by 0°5s. 
If we make this change a permissible smoothing gives the times in Table 5. 

The differences in the times are trivial in relation to the uncertainties 
(+0°3s to 0-5s); but the 1953 times give a considerable change of slope at 15° 
and no other anomaly, while the present ones give a curvature beginning at 14°, 
increasing slightly to 18°-21° and then abruptly decreasing. The observational 
material is slightly in favour of the present solution, and in either case the curvature 
from 14° to 25° is substantially greater than at other distances. 

If Lehmann’s (1934) interpretation of a second arrival (Pd) at about 20° is 
correct, the times should connect with those at 14° by a curve with a downward 
curvature. I take her distances CII for the Azores earthquake (1931 May 20): 
these use the epicentre 37°-6 N, 16°-o W. The earthquake was used in my 
discussion with epicentre 37°-4 N, 15°-9 W, taken from the I.S.S. In a revision 
(not published) I got 37°-44 + 0°-08 N, 16°-17+0°-o5 W. To have a comparison 
as nearly independent of the epicentre as possible I take as trial times of Pd the 
formula 


t = 3 min 19°0s+ 13°70(A— 14°) 


obtained by linear extrapolation from 14°, and for Pr those just found. The theoreti- 
cal interval is then compared directly with the observed one. Since the difference 
of dt/dA for Pr and Pd is much less than either separately, the intervals are much 
less affected by any error of the epicentre than the total times. Lehmann’s 
distances are used. Thus the observed times are earlier than the extrapolated 
ones and permit a small downward curvature for Pd. This was not true for 
comparison with the 1939 table, which made Pr overtake Pd at about 19°-6, 
and made all the Pd observations lie too late to give a downward curvature. The 
stations are in azimuths near 45°, and the corrected epicentre would increase the 
distances by 0°-19+0°-07. The difference of slope would make the calculated 





Table 6 
Observed and calculated differences Pd—Pr 
Pd-Pr 

Station A calc. obs. O-C 

deg s s s 
Strasbourg 20°4 + 8-7 +5°4 —3°3 

Zurich 20°5 + go +4 —5 
De Bilt 20°7 + 96 +5 —4°6 

Chur 20°9 +10°0 +9 -3 
Karlsruhe 210 +10°4 +5°5 —4'9 
Stuttgart ai:3 +11°4 +8 —3°4 
Ravensburg 21°3 +11°4 +8 —3°4 


Mean: —3:9+0°%3 
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difference about o-6s larger in magnitude. If the curvature contributes —4s or 
so at 21° and is due to a term in A, its effects up to 14°. would hardly be detectable. 
But as the amplitudes of Pd are larger than those of Pr, which is itself specially 
large at this distance, it is probable that a considerable curvature for Pd begins 
quite close to 20°. But dt/dA for any part of Pd cannot be less than that for any 
part of Pr, and this from 15° to 16° is already at least 12-8s/deg. It seems that 
the existence of Pd must depend on a rather delicate adjustment of conditions, and 
it would not be at all surprising if it is present for earthquakes in some regions 
and not in others. Lehmann’s reproduced records are convincing evidence that 
there was something to read, whatever the interpretation may be. 

Beyond Ravensburg and up to Copenhagen at 26°-3 Lehmann had 13 stations, 
at none of which was Pd read. It seems that the most probable interpretation is 
that the velocity of P varies continuously with depth; that it first decreases slightly 
to account for the constancy of dt/dA to 14°; and that it then increases a little, 
but in a sufficiently short distance to produce triplication. Then Lehmann’s 
data would refer to the more distant cusp on the time curve. 


2. Japanese earthquakes 


Search was made for earthquakes in the region of Japan with numerous 
observations about 20°, but none was found that was suitable. 


3. North American earthquakes 


The 1939 P table has not been checked for distances under 30° in America; 
the data were for European earthquakes supplemented at large distances by 
observations, mostly in Europe, of earthquakes in other regions. The increase in 
the number of North American stations now makes a check possible. The earth- 
quakes used are as follows. 


1945 Apr. 21, June 30, Oct. 29. 

1946 Feb. 15, Mar. 13 (13 h 21 min and 13h 49 min), May 9, July 18 d 6h. 
1947 Aug. 7d oh. 

1948 Feb. 28. 

1949 Aug. 22, Oct. 31. 


Classification of the readings by quadrants showed few cases where there were 
many observations within 8°. I assume for each quadrant a correction a; and 
corrections to the table 


8° — 12° x1 
12°—15° x2 
15°—17°°5 x3 


20°5—25° x4 
25° —30° X5- 


For each quadrant I made normal equations for a and the five corrections and 
eliminated a. Then the equations for the corrections were combined for all earth- 
quakes and solved. The problem is of course indeterminate, since an addition to 
a and an equal subtraction from all the x s would alter no residual ; and the normal 
equations are essentially for differences of the xs, and at the last stage of the 
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solution there are two equations differing only in -. Linear combinations 
with independent uncertainties are as follows. 


%4 —0°03 X1 —0°27 X2 —0°37 X3 —0°34 %5 = +0°5+0°4 
X2 —O0°29 X1 —0°49 X3— 0°22 X5 = +0°2+0°5 

x3 —0°38 x] —0°62 x5 = +0:0+0°6 

x1 —-%5 = —O°2+1°0 


No corrections are worth making in these ranges of distance. 
For 18° to 20° I get the following means against the other distances. 


deg. x’ d.f. 
18 +0'4 +04 13°8 12 
19 +o°1 + 0°§ a | 10 
20 —o'8 + 0°6 I2°I 5 


(—1-4 +0°7) 3°8 4 


The second row of results for 20° comes from rejection of one entry. The general 
conclusion is that for earthquakes in North America the 1939 table is satisfactory 
to within a second from 8° to 30° and there is no evidence calling for any change. 

For the observations within 8° I fitted corrections a’+ 5A to the residuals 
against the 1953 table. The only useful information for b came from 1946 March 
13, E and W, which gave respectively 


b = —o-15+0-4s/deg, b = —0-05 + 0'9s/deg, 


and in all b = —o-1+0-3s/deg. This is slight evidence that dt/dA at short 
distances agrees better with the European than with the Japanese value. 

Lehmann (1955) has made a study of P and S in N.E. America. The 
observations were not classified by quadrants, and linear forms were fitted up to 
12° by least squares. Thus the uncertainties may be underestimated, and the 
question of a curvature before 12° is undecided. I have treated the observations 
in the same way as the I.S.S. data. There were many more within 8°. Unfortu- 
nately there were not often more than two in the same quadrant for any of the 
ranges of distance used, and a standard error could not be estimated till a general 
solution had been made. 

The observations up to 8° gave standard deviation 1-os on 11 degrees of 
freedom, b = ( +0-15+0-12)s/deg. Combining with the previous data (which 
were given half weight on account of the larger standard deviation) gave 


= (+0711 + 0-11)s/deg. 


This is enough to show that the velocity of P at short distances agrees better with 
the European than with the Japanese value. To reduce as far as possible the effect 
of the uncertainty of b the mean residuals up to 8° were adapted to distance 5°. 

Between 20° and 30° there were very few observations in quadrants where 
they could be usefully compared with other distances. The solutions for x1, x2, x3 
gave the forms with independent uncertainties: 


X1— 0°50 X2—O°50 x3 = +2°5+1°2 
xg— XZ = +1213 


with standard deviation 1-8s on 8 degrees of freedom. 
Because of the much smaller amount of material the uncertainties are much 
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larger than those of the I.S.S. set. It is possible, however, to adopt values from 
the previous set: 


x1 = 0, X2= +02, x3= +01, x4 = +06, x5 = +02, 


and compare the mean residual against these with the mean residual already 
estimated at 5°. Comparison was possible for four quadrants. 


5° weight a+x, weight 


s s 

1935 Nov.1 Ss —o'2 3 —4'9 2 
1939 Oct. 19 WwW +0'9 2 —1'5 2 
1940 Dec. 20 Ww +0°4 4 —63 3 
1952 Apr.9 E +o-2 3 —6°3 6 


The weighted mean difference is —5:48+0-6s; x2 = 4:1 on 3 d.f. This would 
be (within 1s or so) the correction to the 1939 table for distances from 8°-5 to 
30°. The corrections for Europe in this range are between —2s and —3s, so it 
looks as if (relative to 5°) the American times are shorter by about 3s. However, 
Lehmann’s earthquakes were in the eastern half of North America. The I.S.S. 
series were mostly near the west coast, and a possible alternative would be a 
variation within North America. 
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Rock Creep and Thermal Instability 


Harold Jeffreys 


(Received 1958 February 5) 


Summary 


The theory of instability of a fluid heated below can be adapted 
without difficulty to a rock satisfying the modified Lomnitz law of 
creep. It is found that instability is impossible and therefore that 
steady convection currents would not arise. 


The theory of instability of a fluid heated below assumes a finite viscosity, 
and it is desirable to examine whether it needs serious change for rocks obeying 
Lomnitz’s law of creep or my modification (Jeffreys 1958). It turns out that 
the solution for finite viscosity can be adapted with surprising ease. 

In the standard case of viscosity, the temperature disturbance V’ at marginal 
stability satisfies an equation of the form (Jeffreys & Bland 1950): 


G 
VsV’ = —V2V', (1) 
y 


where G depends on various parameters of the problem, including the excess 
temperature gradient over the adiatabic; v, the kinematic viscosity, is made 
explicit. 

Now for a viscous liquid the shear stress is pv times the rate of shear; that is, 
if v is the displacement, the shear stress contains 


) OvK 
SS, 
Prat xi 
For creep, the initial behaviour is similar to that for viscosity, but this stage 
lasts only for a small fraction of a second. At greater times the rate of strain 
decreases, and the stress contains 


w(x + kp) dog Ox, (2) 


with k > 0,0 < a < 1, where p is the Heaviside operator. Hence we get a formal 
adaptation of (1) by replacing v by an operator of the form (u/pp)(1 +kp-*)"1. 
For « = 1 this would express elasticoviscosity; and for general « (1) will be 
replaced by 


V8V' = Hp(1+kp)Vi2V", (3) 


with H = Gp/u. For « = 1 the lowest power of p reproduces (1). 
(1) was found by neglecting all powers of p but the lowest, and hence, since 
only integral powers of p arose, (3) will be correct to the two lowest powers of p. 
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Then the lowest power on the right leads to 


VSV’ = Hkp'veV', - (4) 
and the critical case of a steady motion (p = 0) implies 
VSV’' = o. (5) 


But subject to the boundary conditions of the problem this has no solutions at 
all; for this would be the condition that marginal stability should occur at zero 
temperature gradient. 

Hence under the law of rock creep, heating below cannot produce convection 
currents. 

Physically the reason is as follows. Consider steady convection in the viscous 
case. ‘The excess temperature in the heat columns is maintained against conduc- 
tion by warm liquid flowing in along the bottom, and the velocity is produced by 
shear stresses. ‘The essential point is that constant shear stress in this case implies 
constant rate of shear. But in the case of creep, even if the stresses were constant, 
the rate of shear would tend to zero as time increases ; and hence even if convection 
was started it would not continue to bring in enough heat to replace the loss by 
conduction. 


160 Huntingdon Road, 
Cambridge. 


1958 February. 
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On the Equivalence of Saddle Point Approximations 
and Ray Theory in Elastic Wave Problems 


S. D. Chopra 


(Received 1958 March 24) 


Summary 


It has been known for some time that the Bromwich method of 
expansion in negative powers of exponentials of the formal solutions 
of elastic wave problems for stratified media by wave theory methods 
gives rise to the various reflected and refracted waves predicted by the 
ray theory of geometrical optics. It is shown here by application to an 
example how the saddle point evaluation by the steepest descents method 
of the various terms in such an expansion of the solution leads to exactly 
the same expressions for the displacements in the reflected and refracted 
waves as can be obtained by the geometrical theory. The example 
chosen is that of the elastic waves produced by a harmonic point source 
of dilatational waves situated in the middle of an internal stratum 
bounded on both sides by halfspaces of identical elastic properties, 
all perfectly elastic, homogeneous and isotropic. 

In the course of the demonstration, we give the reflection and 
refraction coefficients for the amplitudes of displacements when a 
plane P or SV wave is incident at an interface between two elastic 
media. These are in the form in which they occur in the wave solutions 
of problems involving reflection and refraction at plane interfaces. 


1. Introduction 


Bromwich (1916, p. 425-428) studied the problem of a bar with one end fixed 
and struck at the other by a particle moving in the direction of the length of the 
bar. He showed that the expansion of the formal operational solution in negative 
powers of exponentials leads to a solution in terms of pulses. The method has 
since been applied by Muskat (1938a, b), Pekeris (1948, p. 46 et seq.), Officer 
(1951), Newlands (1952, p. 222 et seq.), Honda & Nakamura (1953, 1954) and 
others to the formal solutions of elastic wave problems for stratified media. 

Whereas Muskat and others studied problems involving plane boundaries, 
van der Pol & Bremmer (1937) have discussed the propagation of radio waves from 
a spherical point source situated above a spherical earth. They show (1937, 
pp- 844-6; see Bremmer, 1949, Ch. V for further detail) that the saddle point 
approximation for the reflected field gives an expression which may be obtained 
from the primary field by multiplying it by a reflection coefficient corresponding 
to the extra divergence acquired by a thin cone of rays on being reflected from a 
spherical instead of a plane surface. At the same time, the phase factor has to be 
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modified to take account of the distance travelled by the wave. A similar equiva- 
lence of the saddle point approximation with the ray theory has been shown also 
by Gerjuoy (1948, 1953) for a medium consisting of two fluid halfspaces with a 
plane interface and a point source of compressional waves in one of the halfspaces. 

We propose to apply the method used by van der Pol & Bremmer to a problem 
involving reflection and refraction at plane interfaces and show how the saddle 
point evaluation of the various terms in the Bromwich expansion of the wave 
solution is equivalent to the geometric ray-path solution. We give here the dila- 
tational part of the displacement potential 41, in one of the halfspaces, of the distur- 
bance due to a harmonic point source of dilatational waves situated in the middle 
of an internal stratum bounded on both sides by two identical halfspaces with 
which it is in welded contact. As a preliminary to the Bromwich expansion of 4 
in a form which brings out the intimate relationship between the wave solution 
and the geometric theory, we give the reflection and refraction coefficients when 
a plane wave is incident at an interface between two elastic media. We next give 
the expansion of ¢; in terms of these coefficients. Finally, evaluating the displace- 
ments corresponding to two of the terms of this expansion at the saddle point by 
the method of steepest descents and also obtaining the expressions for the cor- 
responding rays by the geometrical theory, we show that the two methods lead to 
identical expressions for each term. In the geometrical derivation of these expres- 
sions we must take account of the fact that when rays are transformed by reflection 
at an interface from the dilatational to the distortional type or vice versa or are 
refracted from one medium into the other, they no longer diverge from a point 
but from a circle round the axis of symmetry. 
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2. Formal solution for a compressional point source in an internal stratum 


We assume the stratum M and the halfspaces M; and Mg to be homogeneous, 
isotropic and perfectly elastic. Let 


p = density, A, u = Lamé’s constants, 
a = {(A+2)/p}* = velocity of dilatational waves, 
B = (u/p)* = velocity of distortional waves, 


for the stratum. The corresponding quantities for the halfspaces M, and Me 
(of identical elastic properties) are denoted by the same symbols with the suffix 1. 
Let 2h be the thickness of the stratum. 

Take the origin at the source O which lies midway between the two interfaces 
and the axis of z perpendicular to the two interfaces with its positive direction 
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pointing into M,. On account of axial symmetry about Oz, we use cylindrical polar 
coordinates r, 0, z. All quantities are independent of 6 so that 0/00 = o. Let 
U, V, W be the displacement components in the 7, 0, z directions. For a point 


source of dilatational waves, V = 0. We introduce displacement potentials ¢, 4 
so that 


be | mb mY 
ae. a ne () 


¢, % determine the dilatational and distortional parts respectively of the disturbance 
and satisfy the equations 


2h 106 ad 1 dd 
—+-—+—=——, (2) 
Or2 or Or -~— Oz? a2 of2 





A, rdop fp Op 1 Ap 
22° - 2 *-* a2” a ae (3) 
r2  roOr r @22 £2 at? 

in M. If 41, 41; $2, %2 denote the corresponding potentials for Mi, Me, then these 
satisfy similar equations with «, 8 replaced by a1, £1 respectively. 

For a harmonic source, the time factor is e”*. This factor will be omitted in 
future from all expressions for potentials and displacements. 

If R = (r? +2), then the potential 49 for a point source of compressional waves 
at O can be represented by Sommerfeld’s integral (Magnus & Oberhettinger, 
1954, P- 34) 


do on ae teRie ie J Jaro “S (4) 


where a = (f2 —w?/«2)# and Jo is the Bessel function of the first kind and zero 
order. The representation (4) is valid for 


—}tn < arga < 3m when —7 < arg(w/x) < o. (5) 


To satisfy the boundary conditions we combine with (4) suitable elementary 
solutions of the differential equations (2), (3) and their analogues for the haifspaces 
M,, Me. ‘These solutions, obtainable by the method of separation of variables, 
are so chosen as to take account of the symmetry about the plane z = o. Then 
the boundary conditions, viz. the continuity of displacements and stresses across 
the interfaces z = +h, lead to four distinct equations. Solving these equations 
we get the values of the unknowns involved in ¢, % etc. Omitting the details of 
these calculations, we obtain 


Ny + Ny! e-2hb) e-ha-z-Wa, d 
re (Ni+ Ni’ e?)e Jor (6) 


, Y1— Yoe-2ha_ Vy e-2hb_ Y,, e-2ha-2ho 





where we use the following notation: 
Ky = wile, Ky = w/B, k, = w/a, Kp = w/ Bi, (7) 


a= (2x2), b= (P—K,P)t, ay = (—K,2)t, by = (—KZ)t, (8) 
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Q = C—hnP, Q = P-hxj, (9) 

Ny = prPa[(uiQs — wl?) + (wQ—pil?)bi), (10) 

Ny! = px Pal(yiQy — wl?)b— (wQ—p0?)bi], (11) 

Yi = p2(Q2 ied C2ab)( C2 = ayb;) _ dupa?’ (aby + ab) - (1 4 

— 2pys€(Q.—ab)(Qy — arby) + y12(C2 — ab)(Q42— Laybr), 

Yo = p%(—Q2—L2ab)({2— arbi) — dupa 423 (abs — ab) — a” 
— 2ppy€2(— Q— ab)(Qy — ayby) + 1% — C2 — ab)\(Q42— Layby) ©? 

Yq = p(Q2 + 2ab)(L2— arbi) — fypan 2x3 (aby — arb) — - 
— 2ppiyl?(O.+ ab)(Q — ayby) + p1%(L2+ ab)(Q42— Payb,) ** 

Ya = p>(Q2—(ab)(l? — aybi) — dupars?y 2(— abi — a1b)— 


- 2p C2(Q _ ab)(Q - ayb) + p1?(C2 - ab)(Q,?2 - {2a,by). (1 5 


The radicals b, a), b; are subject to conditions similar to (5). We shall also use 
the additional notation: 


Ne = px Pla[(u—pr)arb + (u1Qi — pQ)] (16) 
Ng = preg? Cb[(mQ— p21) + (ui — abi] (17) 
Na = pr Pb[(uHQ—prl?)ai + (121 — pl?)a] (18) 
Ns = 2al[p2Q(l? — ayb1) — wpa(Q + £?)(Q4 — aby) + 1*(Q42— aybs)] (19) 


Ne 


2bC[p?Q2(L? — ayby) — pypa(Q + £?)(Q1 — yb) + pr?(Qi?— ayby)). (20) 
The following relations between the Y’s and N’s are easily verified. 
bNs = aNe 
Ni Y3+Ny'Yi = —N3Ns5 | (21) 
Yi Ya+ YoY3 = —Ns5Ne. 


Before proceeding with the expansion of ¢), we give in the next paragraph the 
reflection and refraction coefficients for plane harmonic waves incident at a plane 
boundary between two elastic media and show how the expressions Yj, Yo, Ys, 
Ni, No, Ng, Na, Ns, Ne arise there. The expressions Y4 and Nj’ do not arise directly 
in these coefficients but are connected with the others by means of the relations (21). 


3. Reflection and refraction coefficients for plane waves 


The general equations for reflection and refraction of plane waves at a plane 
interface between two elastic solids were given by Knott (1899, p. 92-96) and Zoep- 
pritz (1919). They did not, however, solve these equations except in a few simple 
cases. Since then many writers have given the results of detailed numerical 
calculations in special cases. Blut (1932), following the method of Zoeppritz, 
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and Muskat & Meres (1940) give the reflection and refraction coefficients for 
plane waves incident at an interface between two elastic media and also give results 
of extensive numerical calculations. The form of their results, however, is not 
suited to our present purpose and it is therefore proposed to obtain these coeffi- 
cients afresh, omitting the details. 

Let H, H, be the two media with their elastic constants denoted by the same 
symbols as for the media M, M, of Section 2. In applying the results of the 
present section, which hold for any pair of media, we identify H, H, with the media 
M, M, respectively. Therefore, here and later, we employ the same symbols as 
in Section 2 to be able to use without change the notation defined in (7)-(20). 


7 














Fic. 2.—Incident P or SV wave and the 
reflected and refracted waves. 


We take the interface between H, Hj as the xOy plane and the z-axis drawn into 
H,. Let motion take place in the xOz plane, and the waves propagate parallel to 
this plane. Then the displacement in the y-direction is zero and all quantities are 
independent of y, so that 0/dy = o. If U, W be the displacement components 
along Ox, Oz respectively, then we can introduce displacement potentials ®, ‘’, 
such that 


o> =F Cl ls 
Ga——p—, Wa —-——-; (22) 
Ox dz dz Ox 
@, give the dilatational and the distortional parts respectively of the motion. 

An incident dilatational (P) or distortional (S) wave in H gives rise, in general, 
to four waves, a reflected dilatational (P’) and a reflected distortional (S’) in H, 
and a refracted dilatational (P,) and a refracted distortional (S;) in Hj. The 
S waves considered here are of the SV type that is, they are polarized in a plane 
perpendicular to the interfaces, but for brevity in writing we omit the V from our 
notation. 

Let e and f denote the angles of incidence and reflection of the P and S waves 
respectively in H, and ¢;, fi, the angles of refraction into Hj. Then the boundary 
conditions—continuity of displacements and stresses across the interface z = o 
—imply Snell’s law 


sine sinf sine, sin fi 


é % 
rae a le (23) 
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w is the circular frequency of the plane harmonic waves under consideration and 
{is defined by (23). We shall identify this f later with the { of Section 2. With ¢ 
given by (23), we also have 


a=ix,cose, b= ixgcosf, a= ik, cose}, 5b = ikg. cosfi, (24) 


where the symbols a, «,, etc. are defined in (7) and (8). Under certain conditions, 
some of the angles e, f, e1, fi; become complex. Later we will also allow £ to become 
complex. In all these cases, the values of the radicals a, b, a, b; remain subject 
to the condition 


—4n < arga, b, aj, by < 4r. (25) 


We can now state our results in terms of the notation of (7)-(20). 

Let Rpp. denote the reflection coefficient for the complex amplitudes of dis- 
placements when P’ is derived from an incident P by reflection. Other reflection 
and refraction coefficients are denoted by similar symbols. Following Byerly 
(1942, p. 164), we adopt the convention that the directions of the displacements 
in the various waves are taken in such a way that all have positive components 
along Ox. These directions are indicated by arrows in Figure 2. 


Incident P wave 


For a P wave incident in the third quadrant in H, suitable values of the poten- 
tials are 


® = A exp{—ix,(x sine+ 2 cose)}+ A’ exp {—ix,(x sine—2zcose)} ) | 
A ; : in H, (26) 
Y = B’ exp{—ix,(x sinf—z cosf)} 
% = A, =pi-t - e1 + 2 cose;)} in Hh. (27) 
¥1 = Bi exp{—ix, (x sin fi + z cosfi)} 


The time factor e“” is omitted in all these. 
Applying the boundary conditions and solving the resulting equations, we 
obtain for the reflection and refraction coefficients 


Yo x Ni 
Rpp. = —, Rpp, = — — 
Yi a Yj 
' (28) 
R ta Ns R ia Ne 
) ee rn, a = ae 
~~ en Ai Yi 
Incident SV wave 
For an SV incident in H, suitable values for the potentials are 
® = A’ exp{—ix,(x sine—2z cose)}, , 
sla re inH, (29) 
Y = Bexp{—ix,(x sinf+z cosf)}+B’ exp {—ix,(x sinf—z cosf)}, 
M,; = A; exp{— ta (* - €1+2 cos é)}, \ in Hh, (30) 
V1 = B, exp{—ix,( x sinfi+2 cosf1)} J 
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and the resulting values of the reflection and refraction coefficients are 





aN ON: 
Rsp. = — Rsp, = = 
aY, val Yi ( ) 
I 
R Ys R san 
ss = —, ss. = ——. 
Yi Bi Yi 


4. Identification of the two definitions of £ 


We now show how the ¢ of (23) may be identified with the { of Section 2 and 
thus prepare for the use of the above reflection and refraction coefficients in the 
Bromwich expansion of 4}. 

Weyl (1919) has shown that, in a halfspace like z > 0, Sommerfeld’s integral 
(4) for the spherical point source may be regarded as the superposition of all plane 
waves of equal intensity whose directions of propagation make with the z-axis 
not only all angles e ly.ug between o and 47 but also a continued series of complex 
angles extending from 47 to 47 +i00. If g be the second polar angle of the direction 
of propagation of any such wave, then the superposition of all these waves gives 
the expression 


im+io 2am 


[exp {—ix,(x sin e cos g+ ysinesing+<2 cos e)} E sine dedg, (32) 
0.47 9 
where E is a constant determining the common intensity of the waves and sin e de dg 
is the element of solid angle. Carrying out the integration with respect to g gives 
im+ioo 
2nE | Jo(rx, sin e) exp (—ix,2 cos e) sine de. (33) 
0,37 
Putting «,sine = { and taking E = «,/2m leads to Sommerfeld’s integral. 
As the substitution «, sine = ( is the same as (23), and both here and in (23), e 
is the angle made by the corresponding plane wave with the z-axis, the ¢ of (23) 
is identified with that of Section 2. 

This identification indicates that we shall also need to consider complex values 
for the angles of incidence, reflection and refraction. In the next paragraph, we 
allow ¢ to assume general complex values. Then e, f, e1, fi will also have general 
complex values, without a simple physical picture. We then define the coefficients 


Rpp. etc. by means of the relations (28) and (31) subject to the condition (25) 
which renders them single-valued functions of C. 


5. Bromwich expansion of ¢ 
The integral in ¢; is of the type 


| CFL) Jo(rh) at, (34) 
0 


where F({) is an even function of ¢. We now regard ¢ as a complex variable and 
Re(w) > 0. We already have Im(w) < o by (5). Making use of the relation 


Jo(re) = tH s0(r£) + dHio(r6), (35) 
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where Hso and H,o are Hankel functions of the first and second kinds respectively 
and of zero order (Jeffreys & Jeffreys 1956, Sec. 21.02), we can transform the path 
of integration from the real axis into the Sommerfeld contour I’ (Figure 3)*, so 
that 


br = [ CR Colro)ae = 4 { CE Hio(ro) a (36) 
0 >A 


I 


[ starts at —ioo in the 4th quadrant and after surrounding the branch lines 
Re(a) = 0, Re(b) = 0, Re(a1) = 0, Re(bi) = © and all poles of F({) to the right of 
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Fic. 3.—Sommerfeld contour I in the {-plane. 
the imaginary axis, returns to —i0o. A typical such pole is shown by a cross in 


Figure 3. For details of the transformation, see Lapwood (1949, p. 73) or New- 
lands (1952, p. 222 et seg.). We thus have 


(Nj he Ny’ e~2hb) e—ha—z—h)a, 


d. 
h=4 [> Hiro. (37) 





Yoe-2ha— Vx e-2hb_ Yi, e-2ha-2hb 


If we take |w| large enough, then |(Yo/Yi)e~2"2|, |(Ys/Yi)e2">| and 
\(Y4/¥i)e-2"2-2hb| are all <1 and we can expand the integrand in ascending 
powers of these exponentials. This expansion is uniformly convergent and we 
can integrate it term by term. We thus get 


* This and other figures are drawn for the case «, >«>{,> , but the argument in the text 
is valid for all cases. 
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eae je7 —ha-z— Wa: Hig(r — arg ™ 





(38) 


making use of the relations (21) and the reflection and refraction coefficients (28), 
(31). This is the required Bromwich expansion of ¢;. We shall denote the succes- 
sive terms in the expansion by 411, ¢12, etc. 

We now show in Sections 6-8 that the evaluation of these terms by the saddle 
point method leads to the same expressions for the displacements as can be ob- 
tained for the corresponding rays by the fully developed geometrical theory. We 
shall illustrate by evaluating by both methods the displacements corresponding to 
11 and $13 only. 


6. Saddle point approximation to the displacement given by ¢1; 
From (38), 


gu = 4 ie — Rees ha~2-has Hio(rl) dl. (39) 


We first replace H;9(r¢) by the first term of its asymptotic expansion, viz. 


2 \t 
Hio(rt) ~ (=) e—tirt—ta) (40) 
are 
valid for —2z < arg (rf) < mand |r| -” (Copson, 1935, p. 336). We then get 
= ( 
du = = sie —Rpp, e~it$-ha—z-ha, dt. (41) 
ed on 
The component Wi; along Oz of the displacement given by 411 is 
; : 
Wes io aoe — 7 OS) E* e-irt—-ha—z—Wa, df, (42) 
z om 


\ 277 oy a 
r 


An account of the method of steepest descents is given by Jeffreys & Jeffreys 
(1956, Section 17.04). For the evaluation of the integral (42), we shall use the 
formula given by them, viz., for large positive r, 
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| X(f) ef dt ~ 
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where Co is the saddle point and « is the angle made with the real axis by the path 
of steepest descents through fo. For its applicability'to (42), we must have r > z. 
We write 





f(%) = —ir$—ha—(z—h)a, (44) 
then the saddle point is given by 
fF, he (z—h)C 
< Sa” @- Kt aa (45) 


Let Co be the root of (45) and eo, e109 the values of e, e; corresponding to { = fo 
as given by (23), then 


Lo Sin €9 SiN €19 
= = (46) 


w a ay 





and (45) becomes 
h tanep+(z—h) taney = r. (47) 


(46) and (47) give real values for the angles eg and ej9 which correspond to the 
minimum-time-path ray PP; of geometrical optics. This ray starts at the source 


c(r,z) 


Zz 
| leon 














Fic. 4.—The Ray PP,. 


O (Figure 4) as a P wave at an angle eg with Oz and after being refracted at A 

at the interface z = h arrives at the point C(r, z) as a P, wave at an angle eyo 

with Oz. Since r >[z, €o Or €19 is close to go° according as a > a or % < a. 
The path of steepest descents through fp is given by the condition that 


a2 

(C— ta)( A is real and negative. (48) 
ot? t=¢, 
We have 

aef i 
(=) = —fha sec? eg + (2z—h)a sec? e;9}. (49) 

ot? t=¢, = 
Let arg w = arg «,, = argk, = —8, so that o < 6 < jn. Then, by condition 


(48), the path of steepest descent through fo makes an angle }7—45 with the 
real axis. The direction of the path of steepest descent through (9 is shown 
by a dotted line in Figure 3. The modification of the path T' to pass through 
in this direction is evidently possible. 
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The saddle point approximation to (42) therefore gives 





eo _ {hseceg (z—A) seceio 
iw sin €9 cos é1o(Rpp )z_¢ exp| — iw + 


a oy 


Wy ~ - 





- (50) 


Since this is the component along Oz of the displacement in a dilatational wave 
which is travelling at C in a direction making an angle ej9 with Oz, the displace- 
ment in the wave is obtained by simply omitting the factor cos ei9 in (50). 


a cos eg[r{h tan eo sec? eg + (z— A) tan e19 sec? e39}]* 











z ‘) 
H a, c(r,z) 

dy c " 

D’ 
A 
G, A’ B Z= h 
dx E 
O, 
dx O,’ 
E’ 

Oo >r 








Fic. 5.—Geometrical derivation of the displacement in the ray PP,. 
LAOG =e, LBOG = e+de, LAO,G, = e,, LBO,G, = e, +de,. 
OH =2z, HC =r =htane+(z—h)tane,, O.,G, =h, OF, =h,+2—h. 
CC’, DD’, AA’, EE’ are arcs of circles about Oz; CD, AE about O,; C’D’, A’E’ about O,’. 


7. Geometrical derivation of the displacement in the ray PP; 


Let a ray starting at the source O (Figure 5) as a P wave at an angle e to Oz 
meet the interface z = h at A, be r fracted into Mj at an angle e; according to 
Snell’s law and be observed at the point C(r, z). Let a neighbouring ray in the 
plane zOA leave the source at an angle e+de to Oz, meet the interface at B, 
and be refracted into Mj at the angle e; +de,. The rays arriving at C and the neigh- 
bouring point D appear to diverge not from O but from O, where the refracted 
rays AC and BD meet. Let O,;G; = hy be the perpendicular from O, to z = A. 
We need the value of 4;. The equation of AC is 


z—h = (r—h tane) cote. (51) 
The position of O; is obtained by solving (51) with the equation 


d 
0 = —hsec®e cotey—(r—h tane) cose? ey -—~ (52) 
e 
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obtained by differentiating (51) with respect to e. We have also 


de, sin e) cose 


Se (53) 


Cos €) sine 


from Snell’s law « sin e; = a sin e. These give 


sin e cos? e; 


b= ———. (54) 


cos? é sin e) 


The P wave starting at O has been derived from the potential 
I . 
bo = = exp(—ix,R), (55) 


where R is the distance from O. Hence the displacement in this wave at A in 
the stratum is given by 


Odo iw (—ie,R) iw RR (56) 
—_—_— NO —— = e ae 9 
OR aR els ah sece i 5 





neglecting the term in 1/R?. This implies that |w| > «/R, or, in other words, the 
distance OA is very large compared to the wave length of the waves under con- 
sideration. We can therefore regard the wave at A as a plane wave and use the 
reflection and refraction coefficients for plane waves given in (28), (31). The 
amplitude of the displacement given by (56) is 


—(iw/ah) cose. (57) 


Since the ray is refracted into Mj, as P;, we have to multiply this by Rpp,, 
so that the amplitude in the ray PP; at A in Mj is 


—(iw/ah) cose. Rpp,. (58) 


On account of axial symmetry about Oz, energy in the medium Mj carried 
along by the ray AC and its neighbouring rays in the PP, wave diverges not from 
the point Oj in a conical element, but in an element of volume formed by revolving 
the area O,CD through a small angle dX about Oz. Hence the amplitudes in the ray 
PP, at A and C are in the limiting ratio of the inverse square roots of the areas 
AEE’ A’ and CDD’C’, since the same amount of energy crosses these areas in unit 
time. We have 


area AEE’A’ = O,Ade,GAdX = h seceyh tane de, dX, (59) 
and 
area CDD'C’ = O\C de; HC dX = (hy + 2—h) sece, rde, dX. (60) 


Hence the amplitude of the displacement at C is approximately 





area AEE’ A’ \* : ; . 
= ( | x amplitude of the displacement at A in Mj 
area CDD’A’ 
iw sine 
=— Rpp,[r{h tane sec? e+ (z—h) tan e; sec? e;}]-+. (61) 


a“ COSé 
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The appropriate phase factor for the displacement at C in the ray PP, is 
. {OA AC _ {Asece (z—h) sece 
exp | -ia( —+—)| = exp | ~i + \|. (62) 


a a4 a ay ° 





Incorporating this into (61) we get the displacement at C in the ray PP}. When 
e, €, are replaced by é, e19 respectively, this gives the same expression as given 


by (50). 


c(r,z) 


Pr 
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O —r 











Fic. 6.—The Ray PSP,. 


8. Displacements in the ray PSP, 


It may be shown as in Section 6 that the displacement corresponding to ¢13 
given by the saddle point evaluation is 





ab . {hsece 2hsecf (z—h)sece 
iw sineRps Rsp, exp | -io| + 2 + \ 
a a1 


, (6 
a cos e[r{h tane sec? e+ 2h tanf sec? f+ (z—h)tan e; sec? e;}]? (93) 





where the angles e, e;, f are now given by the relations 











sine sinf sine 
* - (64) 
ow B O41 
and 
htane+2h tanf+(z—h)tane, = r. (65) 


(64) and (65) give real values for the angles e, f and e; which correspond to the mini- 
mum-time-path ray PSP, of geometrical optics. The path of this ray is shown in 
Figure 6. 

The amplitude of the displacement in the ray PSP, may be derived geometri- 
cally by a repeated application of the method of Section 7. In Figure 7, OABC 
is the path of the ray PSP, and ODEF that of a neighbouring ray in the plane 
2OA. The three-dimensional picture obtained by revolving the figure in the plane 
204A through an angle dX about Oz is not shown in Figure 7 for the sake of clarity. 
The amplitude of the displacement at C is 





iw» (cette eaclael 
ahsece °\O'B.GB) ~~ '\0,C.HC 


iw 








ResRsp 


h’ secf.h tane.h, secey 
ah sece 


R 66 
(h’ + 2h) sec f(hi + z—h) sece;.r (66) 
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But by an immediate extension of (54), 
sin e cos? f sin f cos? e; 


h’ = h, hy = ————-(h' + 2h). 6 
sin f cos*e F sin é) cost - (7) 





Substituting these values in (66), the amplitude of the displacement in the ray 
PSP, is 
iw sine RpsRsp, 





. 68 
a cose[r{h tane sec? e+ 2h tan f sec? f+ (z—h)tane; sec? e;}]* (68) 
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Fic. 7.—Geometrical derivation of the displacement in the ray PSP. 
OH = 2, HC =r = htane+2h tan f+(z—h) tan e. 


The wave PSP, travels the distances OA, AB, BC with the velocities «, 8, «1 
respectively and so the appropriate phase factor in the displacement in this wave 
at C is 





hsece hs z—h) sece 
exp | -io| = x oe _ “\). 


Incorporating this factor into (68), we get the same expression as in (63). 


(69) 


9. Concluding remarks 

The preceding work can be applied also to the other terms in the expansion of 
¢). We thus see that evaluating at the saddle point by the method of steepest 
descents gives the same results as obtained by the fully developed geometrical 
theory. Although we have illustrated the principle of the equivalence of the saddle 
point approximation with the geometrical theory by a single example, the method 
is applicable to all similar problems and the principle is true generally. 

The Sommerfeld integral representation (4) for the field due to a point source 
is obtained by summing up plane waves of equal intensity and suitable directions. 
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The Bromwich expansions of the derived potentials ¢), etc. exhibit these potentials 
as summations of plane waves that can be obtained by all possible reflections and 
refractions of the plane waves contained in (4). Of all the plane waves contained 
in any particular term of the expansion, the saddle point approximation of the term 
picks out that which gives the corresponding minimum-time-path ray of geo- 
metrical optics. We may recall here the remark made by Lapwood (1949, p. 80). 
Evaluating at the saddle point by the method of stationary phase, he gets a result 
which is identical with a result obtained earlier by Jeffreys (1926) for plane waves, 
and observes that this is “..., as it should be, since we have approximated at a 
particular value @ and a particular point ¢.’’ We are here dealing with a fixed 
value of w. 

The importance of the saddle point approximation lies further in the fact that 
at the value of { picked out by it, the phase of the waves forming the integral has 
a stationary value. Here the neighbouring plane waves are nearly in phase and 
add to build up the amplitude. Away from this value of ¢, they fall more and more 
out of phase and there is no building up of the amplitude. The best example of 
this is the expressoin (4) for the spherical point source itself, viz. 


iy [esesaear> ms ao, (70) 
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for which the saddle point approximation gives the exact value 
R- exp(—ix,R). (71) 
The saddle point is given by 


C = x, sine, (72) 


where e is the angle made with Oz by the line joining the source to the point 
(r, z). Thus the value of { at the saddle point gives the direction of propagation 
from the source to the point (r, 2). 

The saddle point approximation to any term in the Bromwich expansion 
gives only a part of the disturbance resulting from that term. To get the complete 
disturbance, we must add the contributions arising from the poles and branch 
lines of the integrand that may have to be passed over in deforming the Sommerfeld 
contour to bring it into coincidence with the path of steepest descent. It is hoped 


to give a fuller account of the problem of the point source in an internal stratum 
at a later date. 
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Notes on Progress in Geophysics 


Some Recent Developments in the Study of Microseisms 
in Great Britain and the United States 


J. Darbyshire and H. M. Iyer 


Summary 


Although microseisms, the small earth tremors of period 4-10s and 
amplitude up to 10m, have been studied for nearly eighty years, many 
of the results and conclusions obtained are still controversial. The 
controversy has been centred mainly round the work done in Britain 
and the United States. An attempt will be made in this article to 
give an account of the present position of the work in both countries 
and to suggest ways and means by which agreement can be reached. 


1. British work 


This work has been largely inspired by the theory put forward by Longuet- 
Higgins (1950). This theory demonstrated how the pressure effect at the sea 
bottom due to a stationary sea wave system could be sufficient to generate micro- 
seisms. Such a system would be formed when trains of waves of the same period 
meet head on and this can happen near the coast when incident and reflected 
waves meet and also near the centre of a fast moving depression where the wind can 
veer round very quickly and trains of waves moving in opposite directions are 
generated. One important consequence of this theory is that the microseism 
period should be half the wave period as was originally observed by Bernard 
(1941). Subsequent British work (Deacon 1947, Darbyshire 1950, 1954) has 
confirmed this ratio but some American results do not fit so well. 

The North Atlantic Ocean with its extensive meteorological coverage is an 
excellent testing ground for the use of microseisms. The results found can then 
be applied to other oceans like the Indian and the Pacific where there are far less 
weather data available. An attempt was thus made to track Atlantic depressions 
using microseisms. Some work on this had already been done by Lee in 1935 
on the basis of microseisms recorded at Kew Observatory. He assumed that the 
microseisms were Rayleigh waves—waves where the ground particles move in 
ellipses contained by the vertical and the direction of propagation. In such a 
simple case, the tangent of the angle of approach of the waves would be given by 
the ratio of the amplitude of the E.W. and N.S. components. The right quadrant 
could be found from the phase relationships between the two horizontal and the 
vertical components. Lee did find some relation between the phase relationships 
and the quadrant of approach. Darbyshire (1954) tried to improve this method 
by assuming that microseisms consisted not only of Rayleigh waves but also of 
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Love waves where the ground particles move in a horizontal plane at right-angles 
to the direction of propagation. It would follow then that the vertical component 
would be free of the effect of Love waves and so records of the vertical and horizon- 
tal components were correlated while they were displaced in time over a definite 
range. From the point of time at which maximum correlation was obtained, the 
phase difference could be determined and from the actual value of the correlation 
coefficients the ratio of Rayleigh to Love wave activity and the angle of approach 
could be estimated. One storm in 1951 was tracked quite successfully by this 
method. Geddes (1957) also tracked a storm with two horizontal seismographs 
by considering the variations in the ratio of the amplitudes, assuming the co- 
existence of Love and Rayleigh waves. 

More recent work has aimed mainly at improving the original correlation 
technique to suit all conditions of storm movement and position (Iyer 1958). 
A new refinement has been the use of band pass filters in analysing microseism 
records so that microseismic activity is divided into different period groups. 
This facilitates the study of microseisms caused by more than one source. A 
detailed’ study using this technique indicates the complexity of the nature of the 
microseismic source. Another interesting result is that the Love waves appear to 
have a random nature like background noise and come from many directions. 
The nature of Love waves and the variation of the Love/Rayleigh ratio with respect 
to the position of the generating source and recording station should hold 
interesting geological information about the nature of the Earth’s crust. 

All this work has not taken into account the effect of refraction of microseisms. 
It is still impossible to allow for this effect when the waves travel over land because 
of lack of knowledge of the elastic properties of the underlying crust but the effect 
while travelling across the sea can be estimated, at least in a qualitative manner, 
if it is assumed that they are waves between the ocean surface and the sea bed of 
the type first suggested by Stoneley (1926). Darbyshire (1956) worked out 
refraction diagrams for Bermuda on this basis and the diagrams explain some of 
the anomalies found by Carder (1955) for this island. Similar diagrams were 
prepared for the British Isles by Darbyshire & Darbyshire (1957) and Iyer & 
others (1958). This work has shown that the effect of refraction could explain 
in many cases the so-called microseism barriers because for certain directions, the 
sea depths form a configuration which acts like a diverging lens on the waves, 
the energy being spread over a wide angle. 

The main development in the instrumental field in Britain has been the 
development by Tucker (1958) of a new horizontal seismograph which has been 
designed primarily for the study of microseisms. It incorporates an electronic 
feed-back circuit so that the instrumental characteristics are determined entirely 
by the constants in the electric circuit and are independent of the mechanical parts. 
The instrument has a flat response from 2 to 10s period, while the magnification 
is 3000. The output can be fed into a pen-recorder. A detailed account of this 
instrument will be published shortly. 


2. American work 


The study of microseisms received a considerable impetus in the United 
States because they could be used to track the hurricanes which advance along the 
Atlantic coast from the Carribbean area. Considerable work on storm tracking 
was done by the Rev. B. Macelwane who used a tri-partite station with horizontal 
Macelwane-Sprengnether seismographs, with accurate timing devices for the three 
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stations. He was very successful in this connection and his methods were extended 
by the U.S. Navy so that between 1944 and 1947 many such tri-partite stations 
were in use. The tri-partite method is not completely satisfactory, however, as 
it is very difficult to identify individual waves on a seismograph trace when the 
distance between stations exceeds 10 miles (Shaw 1921), and for shorter distances, 
measurements of the order of a fraction of a second have to be made on waves of 
about 6s period. A report by Donn & Blaik (1952) suggests that the success of 
the method has not been in proportion to the amount of money and effort spent. 

The indifferent success of the tri-partite method led Marion Gilmore (1953) 
to develop an entirely empirical method. He argued that every storm giving rise 
to microseisms should produce a definite ratio of intensity at two stations depending 
on its own position. He claimed that this principle is independent of the nature, 
origin, composition and path of propagation of the waves. Data collected over 
eight years were used to prepare micro-ratio charts in which iso-ratio lines are 
drawn for pairs of stations. 

The difficulties of this method consist in keeping seismograph characteristics 
constant over a number of years. Moreover, British work has shown that the 
intensity depends not only on the position but also on the strength of the storm 
and its history for some hours previously. 

All the American work so far described has been done by using horizontal 
seismographs and these probably recorded a large proportion of Love wave 
activity which would have the characteristics of noise. 

Future progress in storm tracking would appear to lie in combining the British 
correlation technique with the American tri-partite method. Research in England 
has hitherto been connected mainly with the three-component station working at 
Kew, but another three-component station has now been set up at Cambridge and 
correlation of the vertical and horizontal components of the two stations to ascer- 
tain the phase differences should give promising results. 

The stationary wave theory has never been universally accepted in the U.S.A. 
Carder (1954) and Dinger & Fisher (1955) have produced evidence which supports 
it but other workers, notably, Donn (1954) do not accept the theory and from time 
to time produce evidence against it. This evidence is usually concerned with the 
two to one ratio and Donn concludes that the microseism period is a function of 
the depth of the sea at the storm position. It is, however, very difficult to draw 
any conclusions from observations which do not include both the sea wave and 
microseism spectra. The sea wave spectrum is usually quite wide and there can be 
a selective absorption of wave energy in a narrower frequency range to form micro- 
seisms. There are various ways in which this can be brought about. It follows 
from Longuet-Higgins’s theory that certain frequencies would resonate with various 
depths of water rather like the organ pipe effect. There is also the effect of 
refraction so that waves of one frequency might be concentrated at a particular 
locality. It would be advisable to prepare a refraction diagram for the east coast of 
the United States similar to the one prepared for the British Isles. It is very likely 
that this would explain many of the results found by Donn. Some discrepancies 
from the two to one ratio may also occur due to the action of swell in shallow water 
areas where ordinary progressive sea waves have some effect on the sea bottom. 
This effect would, of course, be of the same period as the swell. This possibility 
was suggested by Bernard (1952). 

American work on the Pacific side has been led by Gutenburg who is also very 
critical of the stationary wave theory. The same objections can be made to this 
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work as in the work of Donn. In the Pacific case, however, the situations are 
even more difficult to deal with as the microseism source is always very far away 
and direct microseisms from the source are very rarely observed. Interesting 
results have been obtained from the three-component station at Palomar which uses 
strain seismographs which can readily distinguish between Rayleigh and Love 
waves. The results confirm the existence of both type of wave in microseisms but 
the Love waves often come from the same direction as the Rayleigh waves and do 
not form a background noise as shown by recent British work. This discrepancy 
may be due to the fact that these microseisms are mainly caused by swell off the 
coast and have not traversed a long ocean path. 

The only theory that can be regarded as a serious rival to the stationary wave 
theory is that of Press & Ewing (1948). This postulates that there are certain 
periods set up between the ocean surface and the sea bed which are associated 
with stationary values of group velocity. These periods could be excited by atmos- 
pheric oscillations of the same period. This system demands, however, a long 
homogeneous path which is very rarely obtained in practice. Although this theory 
would appear quite plausible in some cases, it is not so advanced quantitatively 
as the stationary wave theory. 

Instrumental work in America has been concerned with the development of 
the strain seismographs already mentioned (Benioff 1955). Sharply tuned seismo- 
graphs have also been developed to facilitate the isolation of microseisms from 
different sources ; this is described by Donn, Ewing & Press (1954). 

A description of work carried out in other countries is outside the scope of this 
article but mention might be made of the work of M. Bath in Sweden. He has 
been studying for several years the microseisms recorded at Upsala which he has 
also compared with those at Scoresby Sund and Bergen. He concludes that the 
Upsala microseisms are mainly local in origin, being due to the passage of meteoro- 
logical fronts across the Norwegian coast. He thinks the mechanism of formation 
could be the stationary wave process in some instances and the Press-Ewing 
process in others. The refraction diagrams given by Iyer & others (1958) show 
that no appreciable microseism energy could penetrate into Scandinavia from Atlan- 
tic depressions and this must explain to some extent the very localized effects 
found by Bath. 

Mention has already been made of Bernard who has been using mobile seismic 
stations and taking many records at different places in France. He was the first 
to observe the two to one ratio between waves and microseisms, and his later 
work has confirmed this. Another Frenchman, Imbert, has made a study of micro- 
seisms recorded in Antarctica. 

In conclusion it can be stated that the study of microseisms in Britain is under- 
going a revival in the light of recent developments where full use of modern instru- 
mentation has been made to study the physics of the problem in detail. It is 
hoped that the fresh interest will be shared by other nations and full use made of 
these tiny oscillations as an aid to meteorology, oceanography, and geophysics. 
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Geophysical Discussion of the Royal Astronomical 
Society 


A discussion was held in the rooms of the Royal Astronomical Society on 1958 
January 17 at 16h 15m on Chemical and Radiometric Methods of 
Surveying, under the chairmanship of Sir William Pugh. He remarked on 
the increasing importance of nuclear energy in our economy and on the usefulness 
of radiometric surveying in prospecting for uranium ores, mapping geological 
units, defining structures, and studying mineral deposits. 

The first speaker, Mr F. H. U. Bowie of the Atomic Energy Division of H.M. 
Geological Survey, in introducing his talk on radiometric surveying, pointed out 
the stimulus the atomic bomb has had on prospecting for uranium and thorium 
ores. Before the war these elements were considered only waste products in radium 
extraction. Prospecting is now based on geological knowledge aided by radiation 
counters either of the Geiger-Miiller or scintillation types. Instruments of three 
sorts were described and discussed: portable ones for field use, car-carried record- 
ing meters, and airborne units. Transistorized scintillometers using sodium iodide 
crystals are now the most common and have advantages over the Geiger-Miiller 
tubes of greater statistical accuracy and a smaller volume for the cosmic ray 
background. A portable instrument is chosen on the basis of how well its results 
are reproduceable. The Geological Survey’s Geiger-Miiller tube ratemeter weighs 
3 lb, can detect 5 pulses per second and averages the readings over a short but 
variable time interval. Its reproduceability is +2 per cent. A probe may be used 
with a portable instrument for assaying in the field. An instrument mounted in a car 
may cover a small area but if raised to a height of 54 ft, the area increases to 50 ft 
in diameter. The fall-off of radiation is as distance squared to a height of about 100 
ft but is linear above this due to gamma-ray scatter in the atmosphere. Airborne 
instruments are usually flown at an altitude of 500 ft and cover an area in lane- 
widths of $ mile though for thorough coverage 1,000 ft is preferable. Slides were 
shown giving results of surveys: one compared portable counters over the Dolgelly 
slates and aluvium and a set compared the isorads (line of equal radiation) obtained 
with different flight spacings on an airborne survey of one anomalous area. Bowie 
closed with comments on the geochemical methods of prospecting for uranium 
and thorium ores employed in the Colonies. There, when searching for uranal 
oxides in surface waters, ion exchange resin, colorimetric, or chromatographic, 
analytical methods are used. 

Dr F. S. Webb of the Department of Geology, Imperial College, continued with 
the topic of geochemical prospecting. This aims at locating areas where an ele- 
ment occurs in anomalous amount relative to its average abundance in the Earth’s 
crust. The primary dispersion pattern of such an anomaly relates to the physico- 
chemical processes forming it and includes mineralization and any alterations 
accompanying it. As these may involve the travel of solutions a considerable 
distance along channels, there may be difficulty in locating the source of the 
primary dispersion pattern. A secondary dispersion pattern may result from 
weathering, both chemical and mechanical. Such a dispersion may be residual or 
may be transported by ground and surface waters, soil creep, or biological activity. 
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Dilution presents difficulties for using the analysis of water to detect an anomaly, 
and the study of stream sediment may be more fruitful. (Slides were projected 
showing the Cu, Co and Ni contained in residual soil over a deposit in Uganda, 
and copper analyses of sediment collected along a stream bed.) Plant analyses 
may reflect the plant’s discrimination in absorbing elements from the soil and 
therefore analysis of the soil gives more reliable information. Sensitivity, simplicity 
and speed are the requirements of geochemical methods. Chromatographic and 
colorimetric techniques are the most common, and spectrography and fluorimetry 
have been used in special cases. Geochemical methods are used a great deal for 
mineral prospecting in Canada and Africa. 

The third speaker was Dr 7. E. T. Horne of the Atomic Energy Section of 
H.M. Geological Survey. He spoke on age determination by the products of radio- 
active decay. This may be important in the study of relative age relationships in 
ore genesis. A slide showed how the decay of potassium-40 had been used in 
establishing the epigenetic relationship of an ore from the Colorado Plateau. 
The main emphasis was on uranium-lead-thorium methods. Chemical analysis 
for these three elements plus an isotope study of lead minerals permits the calcula- 
tion of several isotope ratios. Discrepancies may result from loss of material prior 
to the analysis and this be controlled by factors such as grain size or chemical. 
stability. A broad programme of study is often needed to assess these difficulties. 

Dr W. H. Bullerwell of the Geophysics Branch of H.M. Geological Survey 
was the last to speak and his topic was the use of radioactive logging in boreholes. 
A slide comparing one such record with others (self-potential, resistivity, hole 
diameter, and lithology) made from studies in the same borehole, served to intro- 
duce the subject and underline an argument for the correlation of several types of 
logs in making deductions from borehole data. Methods of radioactive logging 
reviewed included those for recording gamma rays and neutrons. The gamma 
ray detection was by the methods Mr Bowie had outlined. A scintillation spectro- 
meter which can give information about the energy of the radiation as well as its 
rate of reception may be useful in determining what elements are its source. 
If a source of neutrons is lowered down a borehole either the secondary neutrons 
or gamma rays, emitted from the walls of the borehole by bombardment from this 
source, may be recorded. A gamma-gamma log may give a measure of the electron 
density which scatters back the gamma-rays from a gamma-source. Development 
of this in conjunction with a scintillation spectrometer is now in progress. 

In the discussion which followed, Dr D. T. Germain-Fones asked Dr Webb 
about the possibility of geochemical prospecting methods being applied to oil 
prospecting. Dr Webb replied that this seemed particularly difficult due to the 
formation of hydrocarbons in the soil from organic matter and not as likely of 
success as in the search for minerals deposit. A question to Dr Webb about the use 
of helicopters in some geochemical prospecting work drew a reply that this was 
for the transport of personnel. Dr V. Armstrong asked Mr Bowie about the cost of 
counters. He said that portable Geiger-Miiller counters cost £45 and scintillo- 
meters about £75. Mr W. A. Willox asked if portable counters were much used 
while making traverses during regional geological surveying. In reply Mr Boyle 
said the use of a scintillometer, weighing only 5 lb by virtue of the use of transistors, 
was being investigated. This might be sensitive enough to contrast geological 
formations but a few feet of overburden could mask the effect. 

The Chairman thanked the speakers and the contributors to the discussion and 
adjourned the meeting at 18h o5m. 
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Book Review 


The Planet Earth 
Edited by D. R. Bates, F.R.S. 
312 pp. 358. net. (1957, Pergamon Press) 


The art of lecturing to, or writing for, an educated audience who know little 
or nothing of one’s own scientific discipline, or of any at all, is usually acquired 
only by hard work and steady application to the waste-paper basket. The most 
successful lecturer or writer begins by explaining the essential facts and principles of 
his subject very simply and directly, avoiding all technical jargon and if he does this 
thoroughly, he will find that subsequently he can carry his audience with him into 
the most advanced ideas. But it requires a great effort of imagination to realize 
how little other people know of facts and ideas that form the half-unconscious 
background of one’s own thinking; and without this effort it will not be possible 
to avoid jargon, to explain each step as it is taken without invoking arguments 
that are familiar to the specialist but to no one else, and to write throughout in 
that simple, straightforward English that demands no effort from the reader but 
every effort from the writer. 

Can there be any doubt, however, that this hard work and hard thinking are 
worthwhile? Apart from the satisfaction of sharing with others one’s own intellec- 
tual interests and the excitement of scientific study, it may be thought that scientists 
have a duty to explain to their patrons, who are now in large part the taxpaying 
public, their aims, methods and results. When, as with the IGY, the work of the 
scientist attracts wide publicity, the duty is imperative. 

The object of this collection of surveys of geophysical studies is stated explicitly 
to be that of providing “the educated public with a background of knowledge 
which will aid them in appreciating the significance of some of the work being 
done”. That being so, a reviewer must principally consider whether the articles 
are, in fact, clear, simple and straightforward, whether each step of an argument 
is explained as it is taken, and whether the author avoids the use of arguments and 
facts which are familiar to the specialist but which only confuse one who is not. 

The reviewer may readily concentrate on this aspect of this book because the 
names of the authors of the articles are themselves guarantees of authority and of a 
thoroughly up-to-date picture. 

In the first chapter, Professor Chapman introduces the IGY. This is followed 
by a discussion of the origin, age and possible ultimate fate of the Earth by Gerard 
P. Kuiper. Professor K. E. Bullen describes the deep interior, Professor J. Tuzo 
Wilson the crust, and Dr G. E. R. Deacon the oceans. The chapter on the geo- 
magnetic field is by E. H. Vestine, while Professor Bates himself contributes that 
on the composition and structure of the atmosphere. Dr E. T. Eady has two 
chapters on climate and on the general circulation of the atmosphere and oceans. 
Dr E. J. Opik writes on ice ages and Dr B. J. Mason on meteorology. Dr J. A. 
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Ratcliffe writes the chapter on the ionosphere, Professor Bates that on the airglow 
and Professor V. C. A. Ferraro one on aurorae and magnetic storms. The chapter 
on meteors is by Professor A. C. B. Lovell, Dr J. G. Wilson writes about cosmic 
rays and Professor J. B. S. Haldane discusses the genesis of life. There is a 
bibliography which contains anything from four to sixteen works for each article, 
some of them popular accounts, others significant papers on the subject. Finally, 
there are six appendices, to explain certain concepts more fully than is appropriate 
in the text. 

Workers in these subjects will find each of the articles well worth reading and 
some of them very valuable containing as they do the latest facts and ideas hot 
from the mint. It is also clear from the list of chapters that almost every 
aspect of the physical constitution and development of the Earth is covered, the 
emphasis being on those subjects most directly involved in the IGY. 

But as essays in popular exposition, these articles are not uniformly successful. 
Some are very good and while assuming little initial knowledge in the reader, 
introduce him to very advanced ideas. The articles of Professor Bullen and 
Dr. Ratcliffe are particularly successful, while Dr Eady’s article on climatology 
gives a good idea of some of the problems of definition and analysis that arise in 
dealing with geophysical data. 

Other chapters, however, are less satisfactory. What, for instance, will the 
intelligent person who is not familiar with celestial mechanics make of the state- 
ment that ““Neptune’s very nearly circular orbit places a million million years as 
the upper limit to the age of the planetary system” in Kuiper’s article or one who 
is not an organic chemist of ‘The polysaccharides consist of hexagonal pyranose 
sugar residues linked by ether linkages in chains which may occasionally branch” 
in Haldane’s? It is probably more than coincidence that these two are much the 
most speculative and are densely packed with ideas and arguments of wide range 
which, however, will be largely unfamiliar to anyone not widely read. Kuiper’s 
chapter may be unfortunately placed at the beginning since it could probably be 
read with more understanding after some of the others. 

The general reader will find some of the articles very good and others not so 
successful. Geophysicists and astronomers will find all the articles interesting 
and some most stimulating. This is a book that should be in all geophysical 
libraries and can be recommended for sixth-form libraries of schools. 


It is well printed and the diagrams are clear and helpful although, for a book of 
this type, not over many. 


A.H.C. 








TEPMUYECKAA MCTOPMA 3EMJIM C Y4YETOM ITEPEMEHHORM 
TEMJIOMNPOBOWHOCTU EE OBOJIOUKRN 


E. A. Jlw6umoea 


JlaHHbIe 0 pacnpocTpaHeHHOCTH paMOaKTMBHBIX JJIEMCHTOB HM TeMJOMNpOBOAHOCTH SemMJM 
TOBOPAT B MOIb3y WpeCTaBeHHA, YTO ee HaYabHaA TemMMepatypa To(r) Oba HuAKe TeM- 
nepaTyphl miaBltenua. B crarbe mpefcTaBleHbl pesyibTaTb pacueTosw T)(r) qua Semin, 
o6pas0BaBmelicA TyTeM OObeMHEHHA YaCTHI, Tas0mbieporo MpoTomMakeTHOrTO OOJIaKa. 
PaccMaTpuBaeTcaA MepepacnpefesteHne pasMOakTMBHBIX DJIEMCHTOB OT OHOPOHOTO K pac- 
COCHHOMY. 

BemjecTBo o60104KM paccCMaTpMBaeTCA Kak JMIJeKTPHK, OONaqaioulMii B OONacTM BBICO- 
KHX TeMMepaTyp He TOJbKO MOJICKYJIAPHOM TenJOMNpOBOAHOCTBIWO, HO TakyKe M JIY4MCTHIM 
Ten000MeHoM. JIA pacueTOB pelleHHA ypaBHeHHA TeMJONPOBOAHOCTH MCHOJNb3yIOTCA 
meToyx dbyHKuun [puna u MeTox ruppaBuMyeckux aHasoruii. 

IlokazaHo, 4TO BO3pacTaHue Ha OonbMOn rpyOuHe oddpeKkTHBHOM TenONpoBOHOCTH, 
o6ycOBIeHHOe TepeHOCoM J1yuNCTON dHeEpruu He NMPMBOAUT K BHAYMTeNbHOM MoTepe BHy- 
TpeHHero Tema. TenuomponoqHOcTb BepxXHux C0eB SemuM OOycOBNeHAa, riaBHbIM OOpa- 
30M, ee MOJIEKYIAPHON CocTaBIAIUeli, KOTOPaA yObIBaeT C pOCTOM TeMMepaTypbI M MoaTOMy 
IpenATCTByeT BHAYMTeIbHOM MoTepe BHYTpeHHero Tena. ITO BefeT K BaKIIOVeHMW, 4TO 
TeJOBYIO UCTOpNIO SeMIM CileyeT pacCMaTpuBaTb Kak BeKOBOe pasorpeBaHue ee Hep 
M OWHOBpeMeHHO CJIab0e OXJIaxKeHMe ee BHEWIHUX CJ10eB. 


CEMCMUYECKHME HABJUOJEHUA B3PbIBOB ATOMHbIX BOMB 
B ABCTPAJIMMU, B 1956 TOY 


B.A. Boam, A.A. Jloiia u JI. JIov. Cammon 
+ + r 


PaccmoTpeubi Ceiicmuyeckne HabJIOqeHuA MpouUcTeKaloumMe OT YeTLIPeX B3PLIBOB Opu- 
TaHCKMX aTOMHEIX 6om6 B Mapaamure (IOnnan Apcrpamnua), B 1956 r. O6e dasnt, Kak P 
Tak u §, perucTpHpoBaMCh Ha MOJeBEIX CTAHIMAX BONb TpanCcaBcTpasMiicKol menesHot 
Joporn, AA KOTOpETX 0,4°< A < 11°, a TaKne Ha Celicmorpade cuctempt Bennosa BOnn3n 
Ajesangbl, IA KOTOporo A = 8°. 

Cpoku nepeyxBuKeHHii yKa3bIBaloT, YTO CKOpocTh Pn = 8,21 + 0,005 KM/cek, a CKOpOCTB 
Sn = 4,75 + 0,01 nm/cek. Mmerrca take yKa3aHuA Ha TO, 4TO CKOpocTH das Pu S, nyTu 
NpOBUsKCHUA KOTOPEIX HAXOMATCA BOMB NOBEPXHOCTH, COOTBeETCTBeHHO OCTUraIT 6,03 + 
0,09 KM/cek uw 3,55 + 0,04 KM/cek. Hukaknx BbInyAKeHHEIX KOJeOaHHii, BOBHMKAaUIMxX 
OnarofqapA MpoxokyeHuio yepes MpoMe*KYTOUHBIe COM, He HaGOao0cb. 

Hadsmio,eHuA MO*KHO OODACHUTL MpeNoOTO+HKMB, YTO CyMleCTByeT BCero OMMH CIO KOPH, 
B KOTOpOM celicmuy¥ecKue CKOpocTH nocToAHHEI. Ilo stoi runotTese TONMMHA KOpHI Nony- 
yaeTcA B 32 KM 110 dase P uw B 39 KM nO chase S. 

B poknaye pasbupariwr1ca Oosee no_poOHnO pesyJbTaTH B3pLIBOB B Mapasmure, paccMoTpeH- 
Hble mpemye Jloiinom. Ha ocHoBaHuM 9TMX aHHBIX Ob OnyOuMKOBaH Hepsi Tpyy, m0 
M3y4eHMW CTPYKTYPbI KOPbI CO CTOPOHEI ABCTpaJIMM, C IPHMeHeHHeM MeTOJa NpeJIOMJIeHHA. 


IMMEPEMEINMEHWE TOBEPXOCTUM OJHOPOJHOTO YITPYrOro 
MOJIYMPOCTPAHCTBA, IITPOMCXO]AMErO BJIATOTAPA TMOT3EMHOMY 
BPAITJAIOWNIEMY MUMITYJIbCY 


X. J. TWerxepuc u HW. M. Jloneman 


VUccnenyerca nepexpuskeHue MOBepXHOCTH OAHOPOAHOTO, yupyroro NoaynpoctrpaHctsBa, 
npoucxoyamero OnarofqapA BpawlawuleMy MMMyJIbCy, TOUKA MpMIO7KeHHMA KOTOPOrO BHYTPH 
ctpyktyph. Ilogpasymepaetca, 4TO OCb Bpallalollero MOMEHTA MapadeibHa NOBepXHOcTH 
M M3MeHeHHe B 3ABMCHMOCTH OT BpeMeHH MOCe_HerO CooTBeTCTByeT yHkuuM H(t) dopmy1 
Xesucaiiza. PesyubTaTh comocTaBuAioTcA c pemenvem IInnnen (1954 r.), pemiapmero Ty 
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me Balla4y MHBIM MeTOROM. OcHOBHOe cBoiicTBO, MpeCTaBIAIOUee MHTepec WIA PusHkKoB, 
oOnapykeHo B OOnacTAX rye HenocpeycTBeHHOM BouHe S npefmectByeT NpeslOMeHHaR 
pouHa SP. Haiigeno, uro B 9THX CIy4AaAX CMeleHMe HA MOBEPXHOCTM HaYMHAeTCA C bee. | 
kxoHeuHol AMIJIMTYAL BO BpeMA NOABIIeCHMA BONHEI SP, a uTo 3a Heli ClefyeT euje Sones 
BbIpavKeHHaA 6eckoneuHocmd Kora NOBepXHOCTH JOCTUTaeT BOTHA S. TakrKe, QA HeOOIB- | 
mux o6sactelt rye BoHa SP He nonBiAeTCA, CMellleHMe HAYMHAeETCA C pe3sKOrO MMMyIbCa 
B Bue dyHKunH § dopmyus Jiupaka. Hpuesse [[mnnes He yKas3biBaloT HM Ha OJHO M3 OTUX 
cBolicts. PesynbTaTH uCCueOBaHMA MOKa3saHbl Ha puc. 2, 3 u 4. 


@MAKTOP BPEMEHIU P JIO 30° 
(Bropaa Cratpsa) 


Tapoavd Jlacediipuc 


VayyaetcaA HOBaA CepHA BeMJeTPACeHH AIA MpoBepku dakTopa ,,P“ qua epponeiicKnx — 
BeMJeTpAceHMi yo 30°, a B ocoGeHHOCTM Mx cocTosHMA OKOO 20°. TlonpapKu Tada | 
apTopa, JjaHHBie B 1939 rogy, xopomlo coBnagawT C NOmpaBKaMM Ha OCHOBaHMM Goes 
paHuMx aHHEIXx 1953 roma. 

HaOmopenua ,,Pd‘“, nmponssezenntie JleiimaH, cooTepecTByioT eé MHTepmperalun, 
Korga ,,Pr“ naronser ,,Pd“ npu npu6nusntTenbHo 14°, HO BO BCAKOM CJly4ae yMBUTeIbHO 
noctonnoctso dt/dA gua ,,Pd“. IlopuqumMomy Bo3MoOxKHO, 4YTO MpOMCXOMT yTpoeHHe nocuep- 
Hero M 4TO HAOMIOAeEHMA OTHOCATCA K OOacTH Gouee OTaNeHHOM TOUKM MepermOa KpuBolt 
BpeMeHH. 

Jina o6nactu Anonum He nOyyeHO NONe3HEIX HOBBIX JaHHbIX. 

HeckoubkO 3emMsleTpACeHu B CeBepHO M WeHTpaubHOH AMepuKe NOATBepKaloT jaH- 
HbIe TaOMIbI 1939 rofa, B palione oT 8° 70 30°, c TOuHOCTBIO JO JOM CeKYH]IbI, 3a HMCKJIN- 
yeHHeM BO3MO?KHOCTH OOaBOuHOK nocTOoAHHOM. JlonouHUTeIbHLIe CBefqeHuA O (hakTope 
BpeMeHH pu paccTOoAHMM MeHee 8°, aHEI B CTaTbe JleiimaH. OTM aHHble yKasbiBaloT, 4TO 
ckopoctTu QA Epponst nour copmagaioT CO CKOPOCTAMM JIA CeBepHoii AMepukKu, HO 4TO 
(bakTOphI BpeMeHM [WIA aMepMKAHCKUX B3eMseTpACeHUit OT 8° Zo 30°, oTHOCHTembHO 5°, OKa- 
3BIBAIOTCA OKONO TBYX CeKYH], KOpoye epponelicKux. 


NOJI3YYECTb TOPHbIX TOPO], U TEPMMUYECKAA HEYCTOMUMBOCTL 
Tapoavd Jlaceddipuc 


Teopuro HeycToHuMBocTu HarpeToi CHHM3Y +KMAKOCTM MO?KHO JIerKO MpPMMeHMTb K Fop- 
HBIM TOpofaM, YOBJICTBOPAIONIMM BUJOMU3SMCHCHHOMY 3aKOHY MOJIsy4ecTu Jlomuura. 

Haiiqeno, 4TO HeyCTOHMYMBOCTL HeEBO3MO7KHA, a CeETOBaTeJIbHO HE BO3SHHMKAIOT KOH- 
BeKIIMOHHbIe MOTOKH. 


O PABHOLWEHHOCTH MPMBJIM;KEHUH, TOJTYYEHHbIX METOJIOM ITEPEBAJIA 
M JIVYEBOM TEOPMU B IIPOBJIEMAX YIIPYTMX BOJIH 


C. J. Yonpa 


Vissectuo, uro Mero Bpomuya, cormacHo KOTOpOMy (opMaJIbHble pelleHua mpoOsem 
yUpyrux BOJIH B CJIOMCTHX Cpeyax, NOJyYeHHbIe MeTOJOM TeOPMM BOJH, pasiaralwTcA B 
CepHlO OTPHMaTeJILHHX BKCHOHEHIMAJIbHBIX YJICHOB, JjaeT OTParKeHHble M MpeJIOMJIeHHHe 
BOJIHLI, WpecCKkasaHHble Jiyyepol TeopHel reoMeTpmuecKoi ONTHKH. 

TlokasaHo Ha MpuMMepe, YTO OLeHKa YICHOB TaKOrO pa3siO*KeHHA MeTOOM nepesama 
yjaeT TOUHIe BLIPAvKeHMA JIA CMeIeHMH OTpasKeHHOM M MpeJOMJIeHHOK BOJIH, NOU yueHHHX 
m3 reoMeTpuyecKoit Teopnn. Vcnoub3s0BaHHEi MpuMep—oToO yupyrne BOJIHEI M3 rapMOoHH- 
ueCKOrO TOYeYHOTO MCTOYUHMKA PaCXOJALIMXCA BOJIH, HaXOfAerocA mocpeyquHe BHYTpeH- 
Hero CJIOAH, PaciOO#KeHHOTO Me?Ky NOJYyMpOCcTpaHCcTBaMM C OMHMMM HM TeMM *Ke YIPyrMMu 
CBOMCTBaMH: UjleasbHOH yIPyrocTbW, OJHOPOHOCTbWO uM Ms0TporHel. 

B pas6ope yaiorcaA TakxkKe KOO PMUMeHTH AMIJIMTY CMeleHui B OTParkeHHOH Mm mpe- 
JOMIIeHHOH BOTHAX, B TOM Ciy4ae, Kora BOHAa mocKkocTn P unum SV nafaeT Ha mepece- 
yeHnMe JByX yUpyrux cpef. OTH KOdPMIMeEHTH B TOM BUe, B KAKOM OHM BCTpPeualwTcA B 
pellleHHAX Bajlau OTpaxKeHHA M NpeJIOMJeHHA BOJIH Ha MepeceyeHuAX MIOCKOCTeN. 
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